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Part 1 
Nonlinear self-adjointness 

1 Preliminaries 

The concept of self-adjointness of nonlinear equations was introduced [1, 3] for 
constructing conservation laws associated with symmetries of differential equa- 
tions. To extend the possibilities of the new method for constructing conservation 
laws the notion of quasi self-adjointness was suggested in [2]. 1 introduce here 
the general concept of nonlinear self-adjointness. It embraces the previous no- 
tions of self-adjoint and quasi self-adjoint equations and includes the linear self- 
adjointness as a particular case. But the set of nonlinearly self-adjoint equations 
is essentially wider and includes, in particular, all linear equations. The construc- 
tion of conservation laws demonstrates a practical significance of the nonlinear 
self-adjointness. Namely, conservation laws can be associated with symmetries for 
all nonlinearly self-adjoint differential equations and systems. In particular, this 
is possible for all linear equations and systems. 

1.1 Notation 

We will use the following notation. The independent variables are denoted by 

X = {x^, . . . 

The dependent variables are 

u^{u\...,u'^). 

They are used together with their first-order partial derivatives : 

^(i) = {<}, <-AK), 

and higher-order derivatives U(^2), ■ ■ ■ , ^^(s), ■ ■ ■ , where 

«w = K...J, <...,^ = A,---A.K). 

Here Dj is the total differentiation with respect to a;* : 

A locally analytic function /(x, -u, . . . of any finite number of the 

variables x, u, ^(2), ... is called a differential function. The set of all differen- 
tial functions is denoted by A. For more details see [4], Chapter 8. 
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1.2 Linear self-adjointness 



Recall that the adjoint operator F* to a linear operator F in a Hilbert space H 
with a scalar product {u, v) is defined by 

{Fu, v) = {u, F*v) , u,v eH. (1.2) 

Let us consider, for the sake of simplicity, the case of one dependent variable 
u and denote by H the Hilbert space of real valued functions u{x) such that u^{x) 
is integrable. The scalar product is given by 



{u,v) = / u{x)v{x)dx. 

Let F be a linear differential operator in H. Its action on the dependent variable 
u is denoted by F[u\. The definition (1.2) of the adjoint operator F* to F, 

{F[u],v) = {u,F*[v]), 

can be written, using the divergence theorem, in the simple form 

vF[u] - uF*[v] = Di{p'), (1.3) 

where visa new dependent variable, and are any functions of x, u, v, , , . . 
It is manifest from Eq. (1.3) that the operators F and F* are mutually adjoint, 

{F*y = F. (1.4) 

In other words, the adjointness of linear operators is a symmetric relation. 

The linear operator F is said to be sclf-adjoint if F* = F. In this case we say 
that the equation F[u] = is sclf-adjoint. Thus, the self-adjointness of a linear 
equation F[u] =0 can be expressed by the equation 

F*[v] =F[u]. (1.5) 

v=u 

1.3 Adjoint equations to nonlinear differential equations 

Let us consider a system of m differential equations (linear or nonlinear) 

Fa{x,u,U{i),...,U{s)) ^0, Q; = l,...,m, (1.6) 

with m dependent variables u = {u^, . . . ,u"^). Eqs. (1.6) involve the partial 
derivatives . . . , U(^s) up to order s. 



Definition 1.1. The adjoint equations to Eqs. (1.6) are given by 

F*{x,u,v,U(i),V(i),...,U(s),v^s)) =(^, a = l,...,m, (1.7) 

with 

F*{x,u,v,u^i),v^i), . . . ,u^s),v^s)) = j^, (1-8) 
where C is the formal Lagrangian for Eqs. (1.6) defined by ^ 

m 
/9=1 

Here v — {v^, . . . , v'^) are new dependent variables, . . . , V{^s) are their deriva- 
tives, e.g. i'(i) = {^i}, '^i — Di{v°'). We use 6/Su°' for the Euler- Lagrange 
operator 

S c) d 



du^" ' du'^ , 

S=l 



SO that 



5{vPF,) _ djv^F,) ( d{v^F,) \ fdiv^F, 



The total differentiation (1.1) is extended to the new dependent variables: 

The adjointness of nonlinear equations is not a symmetric relation. In other 
words, nonlinear equations, unlike the linear ones, do not obey the condition (1.4) 
of mutual adjointness. Instead, the following equation holds: 

(F*)* = F (1.11) 

where F is the linear approximation to F defined as follows. We use the tem- 
porary notation F[u] for the left-hand side of Eq. (1.6) and consider F[u + w] 
by letting w <^ 1. Then neglecting the nonlinear terms in w we define F by the 
equation 

F[u + w] ^ F[u] + F[w] (1.12) 

For linear equations we have F = F, and hence Eq. (1.11) is identical with Eq. 
(1.4). 



^See [1]. An approach in terms of variational principles is developed in [5]. 
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Let us illustrate Eq. (1.11) by the equation 

F = Uxy — smu — 0. (1-13) 

Eq. (1.8) yields 

F* = -^[v{uxy — sinti)] = — vcosu (1-14) 

and 

(F*)* = -^[w{vxy — vcosu)] — Wxy — wcosu. (1-15) 

Let us find F by using Eq. (1.12). Since sinw w, cosw 1 when w <^ 1, we 
have 

F[u + w] = {u + w)xy — sin(-u + w) 

= Uxy + Wxy — sin u COS w — sin w cos u 
~ Uxy — smu + Wxy — w COS u, 

— F[u] + Wxy — W COS u. 

Hence, by (1.12) and (1.15), we have 

F[w] ^ Wxy - wcosu ^ (F*)* (1.16) 
in accordance with Eq. (1.11). 

1.4 The case of one dependent variable 

Let us consider the differential equation 

F(x,ii,ii(i), . . . = (1.17) 

with one dependent variable u and any number of independent variables. In this 
case Definition 1.1 of the adjoint equation is written 

F*{x,u,v,u^i),v^i),...,u^s),V{s)) =0, (1.18) 

where 

F*{x,u,v,U(^i),v^i), . . . ,u^s),V(s)) = • (1-19) 



1.5 Construction of adjoint equations to linear equations 

The following statement has been formulated in [1, 3]. 

Proposition 1.1. In the case of linear differential equations and systems, the 
adjoint equations determined by Eq. (1.8) and by Eq. (1.3) coincide. 

Proof. The proof is based on the statement (see Proposition 7.1 in Section 7.2) 
that a function Q{u,v) is a divergence, i.e. Q = Di{h^), if and only if 



5Q_ 



= 0, 



6Q_ 
5v°' 



— 0, a — 1, ... ,171. 



(1.20) 



Let the adjoint operator F* be constructed according to Eq. (1.3). Let us consider 
the case of many dependent variables and write Eq. (1.3) as follows: 

v^F[s[u]=u^F;[v] + D,{p'). (1.21) 

Applying to (1-21) the variational differentiations and using Eqs. (1-20) we obtain 

Hence, (1.8) coincides with F*[v] given by (1.3). 
Conversely, let F*[v] be given by (1.8), 



Consider the expression Q defined by 



S{v^F^[u]) 



Q = v^F,[u] - u^F;[v] ^ v^F,[u] - • 

Applying to the first expression for Q the variational differentiations S/Su"" we 
obtain 

Applying 5/5v°' to the second expression for Q we obtain 
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The reckoning shows that 



u 



Su^ 



Fa[u] 



Sv°' 



Su^ 



u 



A^'F,[u]) 
5u^ 



^F^[u]. 



(1.22) 



Thus Q solves Eq. (1.20) and hence Eq. (1.21) is satisfied. This completes the 
proof. 
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Remcirk 1.1. Let us discuss the proof of Eq. (1.22) in the case of a econd-order 
hnear operator for one dependent variable: 

F[u\ — a^^{x)uij + U'{x)ui + c{x)u. 

Then we have: 

5{yF[u]) 



u- 



^u[cv- vDi{b') + vDiDj{a'^) - Uvi + 2viDj{a'^) + a'^Vij] . 



Su 

Whence, after simple calculations we obtain 

S 



Sv 



S{vF[u]y 
Su 



[cu + b'ui + a'^Uij] + {DiDj{a'^u) - Di{a'^Uj) - Di[uDj{a'^)]} 



and, noting that the expression in the braces vanishes, arrive at Eq. (1.22). 

Let us illustrate Proposition 1.1 by the following simple example. 
Example 1.1. Consider the heat equation 

F[u] =ut-u^^ = (1.23) 
and construct the adjoint operator to the linear operator 

F^Dt-Dl (1.24) 

by using Eq. (1.3). Noting that 

vut = Dt{uv) - uvt, 

VUxx = Dx{vUx) — VxUx = Dx{vUx — UVx) + UVxx 

we have: 

vF[u] = v{ut - Uxx) = u{-vt - Vxx) + Dt{uv) + Dx{uvx - VUx). 

Hence, 

vF[u] - u{-vt - Vxx) = Dt{uv) + Dx{uvx - VUx). 

Therefore, denoting t = a;^, a; = x^, we obtain Eq. (1.3) with F*[v] = —Vt — Vxx 
and = uv, = uVx — vUx- Thus, the adjoint operator to the linear operator 
(1.24) is 

F* = -Dt - Dl (1.25) 
and the adjoint equation to the heat equation (1.23) is written —vt — Vxx = 0, or 

vt + vxx = ^. (1.26) 

The derivation of the adjoint equation (1.26) and the adjoint operator (1.25) 
by the definition (1.19) is much simpler. Indeed, we have: 

5{vut-vuxx) ^, . ^.2/ N / , ^ 
F = = -Dt[v) - D^[v) ^-[vt + Vxx)- 
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1.6 Self-adjointness and quasi self-adjointness 

Recall that a linear differential operator F is called a self-adjoint operator if it is 
identical with its adjoint operator, F — F*. Then the equation F[u\ = is also 
said to be self-adjoint. Thus, the self-adjointness of a linear differential equation 
F[u\ = means that the adjoint equation F*[v\ = coincides with F[u] = upon 
the substitution v = u. This property has been extended to nonlinear equations 
in [1]. It will be called here the strict self-adjointness and defined as follows. 

Definition 1.2. We say that the differential equation (1.17) is strictly self- adjoint 
if the adjoint equation (1.18) becomes equivalent to the original equation (1.17) 
upon the substitution 

v^u. (1.27) 

It means that the equation 

F*{x,U,U, . . . = \F{x,U, . . . (1.28) 

holds with a certain (in general, variable) coefficient A. 
Example 1.2. The Korteweg-de Vries (KdV) equation 

is strictly self-adjoint [3]. Indeed, its adjoint equation (1.18) has the form 

Vt = Vxxx + "^Vx 

and coincides with the KdV equation upon the substitution (1.27). 

In the case of linear equations the strict self-adjointness is identical with the 
usual self-adjointness of linear equations. 

Example 1.3. Consider the hnear equation 

Utt + a{x)uxx + b{x)ux + c{x)u = 0. (1-29) 
According to Eqs. (1.18)-(1.19), the adjoint equation to Eq. (1.29) is written 

-^{v[utt + a{x)uxx + b{x)ux + c{x)u\) = D^{v) + Dl{av) — Dx{bv) -\- cv — 0. 

Upon substituting v = u and performing the differentiations it becomes 

Utt + auxx + (2a' - b)ux + (a" - b' + c)u = 0. (1.30) 

According to Definition 1.2, Eq. (1.29) is strictly self-adjoint if Eq. (1.30) coin- 
cides with Eq. (1.29). This is possible if 

h{x)^a\x). (1.31) 
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Definition 1.2 is too restrictive. Moreover, it is inconvenient in the case of 
systems with several dependent variables u = (u^, . . . ,u"^) because in this case 
Eq. (1.27) is not uniquely determined as it is clear from the following example. 

Example 1.4. Let us consider the system of two equations 

1,2 2 2 A 

Uy + u u^-Ut ^0, 

ul-ul = (1.32) 

with two dependent variables, u = [u^, u'^), and three independent variables t, x, y. 
Using the formal Lagrangian (1.9) 

C = v\ul + u\l-u^,) + v\ul-ul) 

and Eqs. (1.8) we write the adjoint equations (1.7), changing their sign, in the 
form 

vl + u\\ - v\ = 0, 

v\-vl^^. (1.33) 
If we use here the substitution (1.27), v = u with v = (v^, v^), i.e. let 

— u^, — u'^, 
then the adjoint system (1.33) becomes 

ul + u\l - ul = 0, 

which is not connected with the system (1.32) by the equivalence relation (1.28). 
But if we set 

y-^ = u^, f ^ = u^, 
the adjoint system (1.33) coincides with the original system (1.32). 

The concept of quasi self-adjointness generalizes Definition 1.2 and is more 
convenient for dealing with systems (1.6). This concept was formulated in [2] as 
follows. 

The system (1.6) is quasi self-adjoint if the adjoint system (1.7) becomes 
equivalent to the original system (1.6) upon a substitution 

V = cp{u) (1.34) 

such that its derivative does not vanish in a certain domain of u, 

cp'{u)^0, where cp' (u) = \\^^^\\. (1.35) 
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Remcirk 1.2. The substitution (1.34) defines a mapping 

v"" — ip°'{u), Q; = l,...,m, 

from the m-dimensional space of variables u — {u^, . . . , u'^) into the m-dimensional 

space of variables v = {v^, . . . ,v"^). It is assumed that this mapping is contin- 
uously differentiable. The condition (1.35) guarantees that it is invertible, and 
hence Eqs. (1.7) and (1.6) are equivalent. The equivalence means that the fol- 
lowing equations hold with certain coefficients : 

F*{^^'^^'P^---^'^{s),'P{s)) ^ ^iFp{x,u,...,u^s)), a^l,...,m, (1.36) 
where 

ip^{ip-{u)}, (^(,) = {Ai---A.(¥'"(^i))}, a^l,...,s. (1.37) 
It can be shown that the matrix ||A^|| is invertible due to the condition (1.35). 
Example 1.5. The quasi self-adjointness of nonhnear wave equations of the form 

Utt - Uxx = fit, X, U, Ut, Uro) 

is investigated in [6]. The results of the paper [6] show that, e.g. the equation 

Utt - Uxx + - ul = (1.38) 

is quasi self-adjoint and that in this case the substitution (1.34) has the form 

V = e". (1.39) 

Indeed, the adjoint equation to Eq. (1.38) is written 

vtt - Vxx - '^vutt - 2utVt + 2vuxx + ^UxVx = 0. (1-40) 

After the substitution (1.39) the left-hand side of Eq. (1.40) takes the form (1.36): 

Vtt - Vxx - '^vutt - 2utVt + 2vuxx + 2uxVx = -e"[titt - Uxx + - ul]. (1.41) 

It is manifest from Eq. (1.41) that v given by (1.39) solves the adjoint equation 
(1.40) if one replaces u by any solution of Eq. (1.38). 

In constructing conservation laws one can relax the condition (1.35). There- 
fore I generalize the previous definition of quasi self-adjointness as follows. 
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Definition 1.3. The system (1.6) is said to be quasi self-adjoint if the adjoint 
equations (1.7) are satisfied for aU solutions u of the original system (1.6) upon 
a substitution 

v''^(p''{u), a = l,...,m, (1.42) 

such that 

ifiu) ^ 0. (1.43) 

In other words, the equations (1.36) hold after the substitution (1.42), where not 
all (fi°'{u) vanish simultaneously. 

Renicirk 1.3. The condition (1.43), unlike (1.35), does not guarantee the equiv- 
alence of Eqs. (1.7) and (1.6) because the matrix ||A^|| may be singular. 

Example 1.6. It is well known that the linear heat equation (1.23) is not self- 
adjoint {not strictly self-adjoint in the sense of Definition 1.2). It is clear from 
Eqs. (1.23) and (1.26). Let us test Eq. (1.23) for quasi self-adjointness. Letting 
V = ip{u), we obtain 

vt = (p'ut, = (p'uj;, v^^ = (p'u^^ + (p"ul, 

and the condition (1.36) is written: 

(p'{u)[ut + u^j] + (p"{u)ul ^X[ut- U^x]- 

Whence, comparing the coefficients of Ut in both sides, we obtain A = ^'{u). Then 
the above equation becomes 

^'{u)[ut Uxx] + <P"iu)ul = (p'{u)[ut - U^x]- 

This equation yields that (p'{u) — 0. Hence, Eq. (1.23) is quasi self-adjoint with 
the substitution v — C, where C is any non- vanishing constant. This substitution 
does not satisfy the condition (1.35). 

Example 1.7. Let us consider the Fornberg-Whitham equation [7] 

Ut - utxx - uu^xx - 'iu^u^.x + uux + Ux^Q. (1-44) 

Eqs. (1.18)-(1.19) give the following adjoint equation: 

F* = -Vt + vtxx + uvxxx - uvx -Vx = 0. (1-45) 

It is manifest from the equations (1.44) and (1.45) that the Fornberg-Whitham 
equation is not strictly self-adjoint. Let us test it for quasi self-adjointness. In- 
serting in (1.45) the substitution v = <f{u) and its derivatives 

Vt = (f'ut, Vx = (p'Ua:, Va:x = v'Uxx + v"ul, Vtx = v'utx + ^"utU^, ... , 

then writing the condition (1.36) and comparing the coefficients for Ut, Utx, u^x, ■ ■ ■ 
one can verify that (p'{u) = 0. Hence, Eq. (1.44) is quasi self-adjoint but does not 
satisfy the condition (1.35). 
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2 Strict self-adjointness via multipliers 



It is commonly known that numerous linear equations used in practice, e.g. lin- 
ear evolution equations, are not self- adjoint in the classical meaning of the self- 
adjointness. Likewise, useful nonlinear equations such as the nonlinear heat equa- 
tion, the Burgers equation, etc. are not strictly self-adjoint. We will see here that 
these and many other equations can be rewritten in a strictly self- adjoint equiv- 
alent form by using multipliers. The general discussion of this approach will be 
given in Section 3.7. 

2.1 Motivating examples 

Example 2.1. Let us consider the following second-order nonlinear equation 

Uxx + f{u)ux -Ut^O. (2.1) 

Its adjoint equation (1.18) is written 

Vxx-f{u)v, + vt = 0. (2.2) 

It is manifest that the substitution v = u does not map Eq. (2.2) into Eq. (2.1). 
Hence Eq. (2.1) is not strictly self-adjoint. 

Let us clarify if Eq. (2.1) can be written in an equivalent form 

lj{u) [u^^ + f{u)u^ -ut\^Q (2.3) 

with a certain multiplier ii{u) 7^ so that Eq. (2.3) is strictly self-adjoint. The 
formal Lagrangian for Eq. (2.3) is 



C = viJ,{u)[u:cx + fiu)u^ - Ut\. 



We have: 



5jd 

— = Dl[fi{u)v] - aMu)fiu)v] + Dt[iiiti)v] 
+ ii'{u)v[u^^ + f{u)u^ - ut] + ii{u)f'{u)vu^ , 

whence, upon performing the differentiations, 

SjC 
ou 



The strict self-adjointness requires that 



6C 
5u 
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This provides the following equation for the unknown multiplier ii{u) : 

{jji + 2ufi')u^^ + {2fi' + uii")ul - fifu^ + fiUt = X[u^x + f{u)ux - Ut]- (2.4) 

Since the right side of Eq. (2.4) does not contain ul we should have 2n' + un" = 0, 
whence /i — Ciu"^ + Cg. Furthermore, comparing the coefficients of Ut in both 
sides of Eq. (2.4) we obtain A = — Now Eq. (2.4) takes the form 

(C2 - Ciu~^)u^^ - {Ciu~^ + C2)fu^ = -(CiM"^ + C2)[u,^x + /(m)mJ 

and yields C2 = 0. Thus, = Ciu^^. We can let Ci = — 1 and formulate the 
result. 

Proposition 2.1. Eq. (2.1) becomes strictly self-adjoint if we rewrite it in the 
form ^ 

- [ut - u^x - f{u)u^] = 0. (2.5) 



Example 2.2. One can verify that the nth-order nonlinear evolution equation 

(2.6) 



dt ^ ' dx"" J \ )^ ^ 



with one spatial variable x is not strictly self-adjoint. The following statement 
shows that it becomes strictly self-adjoint after using an appropriate multiplier. 

Proposition 2.2. Eq.(2.6) becomes strictly self-adjoint upon rewriting it in the 
following equivalent form: 

1 



uf{u) 



0. 



(2.7) 



Proof. Multiplying Eq. (2.6) by /i('u) and taking the formal Lagrangian 

C = Vn{u)[Ut - f{u)Un], 

where Un — D2{u), we have: 

— = -Dt[(x{u)v] - D!^[fi{u)f{u)v] + vn'{u)ut - v[/x(m)/(m)]'m„. 

Noting that —Dt[ji{u)v] + v^'{u)ut — —^{u)vt and letting v — uwe obtain 
5C 



Su 



= -H{u)ut - D^[id{u)f{u)u] - [id{u)f{u)]'uun. 



If we take ii{u) — [uf{u)] ^, then ii{u)f{u)u — 1, ii{u)f{u) — u ^, and hence 



Su 



uf(u) 



[ut - f{u)Un]. 



Thus, Eq. (2.7) satisfies the strict self-adjointness condition (1.28) with A = — 1. 
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2.2 Linear heat equation 

Taking in (2.5) f{u) = 0, we rewrite the classical linear heat equation ut — u^x 
in the following strictly self-adjoint form: 

- k - Uxx] = 0. (2.8) 
u 

This result can be extended to the heat equation 

Ut-Au = 0, (2.9) 

where Am is the Laplacian with n variables x — {x^, . . . , x"). Namely, the strictly 
self-adjoint form of Eq. (2.9) is 

-[ut- Au]=0. (2.10) 
Indeed, the formal Lagrangian (1.9) for Eq. (2.10) has the form 

V 

C = — [Ut — Am . 
u 

Substituting it in (1.19) we obtain 

f- = -a(-)-a(^)-4i«,-a«]. 

Upon letting v — u it becomes 

F*^--[ut-Au]. 
u 

Hence, Eq. (2.10) satisfies the condition (1.28) with A = —1. 

2.3 Nonlinear heat equation 

Consider the nonlinear heat equation Ut — {k{u)ux) = 0, or 

Uf — k{u)uxx — k'{u)ul — 0. 

Its adjoint equation has the form 

Vt + k{u)vxx = 0. 

Therefore it is obvious that (2.11) does not satisfy Definition 1.2 
strictly self- adjoint if we rewrite it in the form 

- [ut - k{u)uxx - k'{u)u^^ = 0. 
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(2.11) 



. But it becomes 
(2.12) 



Indeed, the formal Lagrangian (1.9) for Eq. (2.12) is written 

V 

jC^ - [ut- k{u)u^x - k'{u)ul] . 
Substituting it in (1.19) we obtain 

F* = - A - Dl k{u)) + 2A k'{u)u,) 

- - k'{u)Uxx - - k"{u)ul - ^ [ut- k{u)u^^ - k'{u)ul] . 
Letting here v — uwe have: 



F* = -- [ut - k{u)uxx - k'{u)ul] . 
u 

Hence, Eq. (2.10) satisfies the strict self-adjointness condition (1.28) with A = — 1. 

2.4 The Burgers equation 

Taking in (2.5) f{u) = u we obtain the strictly self-adjoint form 

- [ut - Uxx] -u^^O (2.13) 

u 

of the Burgers equation ut = Uxx + uux- 

2.5 Heat conduction in solid hydrogen 

According to [8], the heat conduction in solid crystalline molecular hydrogen 
at low pressures is governed by the nonlinear equation (up-to positive constant 
coefficient) 

Ut = u^Au. (2.14) 
It is derived from the Fourier equation 

dT 

pc,-^ = V -{kyT) 

using the empirical information that the density p at low pressures has a con- 
stant value, whereas the specific heat and the thermal conductivity k have the 
estimations 



It is also shown in [8] that the one-dimensional equation (2.14), 

Ut^u^Uxx, (2.15) 
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is related to the linear heat equation by a non-point transformation (Eq. (5) in 
[8]). A similar relation was found in [9] for another representation of Eq. (2.15). 
The non-point transformation of Eq. (2.15) to the linear heat equation 

Ws = w^^ (2.16) 

is written in [10] as the differential substitution 

t — s, X = w, u = w^. (2-17) 

It is also demonstrated in [10], Section 20, that Eq. (2.15) is the unique equation 
with nontrivial Lie-Backlund symmetries among the equations of the form 

ut^ f{u) + h{u,Ux), f'{u)y^O. 

The connection between Eq. (2.15) and the heat equation is treated in [11] as a 
reciprocal transformation [11]. It is shown in [12] that this connection, together 
with its extensions, allows the analytic solution of certain moving boundary prob- 
lems in nonlinear heat conduction. 

Our Example 2.2 from Section 2.1 reveals one more remarkable property of 
Eq. (2.15). Namely, taking n ~ 2 and f{u) ~ in Eq. (2.7) we see that Eq. 
(2.15) becomes strictly self-adjoint if we rewrite it in the form 

^ = ^. (2.18) 
u-^ u 

2.6 Harry Dym equation 

Taking in Example 2.2 from Section 2.1 n = 3 and f{u) — we see that the 
Harry Dym equation 

Ut — u^Uxxx = (2-19) 
becomes strictly elf-adjoint upon rewriting it in the form 



U 



= 0. 



2.7 Kompaneets equation 

The equations considered in Sections 2.1 - 2.6 are quasi self-adjoint. For example, 
for Eq. (2.6) we have 

F* = -vt-Dl{f{u)v)-vf\u)un, 
whence making the substitution 

1 

V — 
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we obtain 

r.* r f f 



^Ut- —Un^ -^[Ut- f{u) Un]. 



P f f 

Hence, Eq. (2.6) is quasi self-adjoint. 
Example 2.3. The Kompaneets equation 



ut^ ^D,[x\u^ + u + u^)] (2.20) 



provides an example of an equation that is not quasi self-adjoint. Indeed, Eq. 
(2.20) has the formal Lagrangian 

C — v[—Ut + x^Uxx + {x^ -\- Ax -\- 2x^u)ux + 4:x{u + u^)]. 

The calculation yields the following adjoint equation to (2.20): 

SjO 

— = vt + x^v^x - x^{\ 2u)vx 2{x 2xu - \)v = 0. (2.21) 



Letting v — (p{u) one obtains: 

= '{)\u)\ut x^u^x - x^(l 2u)u^ 



6C 

6u 



v=(p{u) 

2„,2 



+ (p"{u)x^ul + 2{x + 2xu - l)ip{u) 
Writing the quasi self-adjointness condition (1.36) in the form 
5C 



5u 



— \[—Ut + X Uxx + {x + 4x + 2x u)ux + 4:x{u + u )] 

v=(f{u) 



and comparing the coefficients for Ut in both sides one obtains A = —(p'{u), so 
that the quasi self-adjointness condition takes the form 

(p'{u)[ut + x^Uxx — x^{^ + 2u)u^ -|- (fi"{u)x^ul + 2{x + 2xu — l)ip{u) 
= (f'{u)[ut — x^Uxx — {x^ -|- 4a; -|- 2x^u)Ux — 4a;(-u + u^)]. 

Comparing the coefficients for u,ji.,ji. in both sides we obtain ^'{u) = 0. Then the 
above equation becomes {x + 2xu — l)(p{u) = and yields (p{u) = 0. Hence the 
Kompaneets equation is not quasi self-adjoint because the condition (1.43) is not 
satisfied. 

But we can rewrite Eq. (2.20) in the strictly self-adjoint form by using a more 
general multiplier than above, namely, the multiplier 

x'^ 

pL^ — ■ (2.22) 
u 
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Indeed, upon multiplying by this /i Eq. (2.20) is written 

— Ut^-D^ [x^{ux + U + u^)] . 
u u ^ 



Its formal Lagrangian 

u 



satisfies the strict self-adjointness condition (1.28) with A = — 1 



5_C 

Su 



u 



Remark 2.1. Note that v = x"^ solves Eq. (2.21) for any u. The connection 
of this solution with the multiplier (2.22) is discussed in Section 3.7. See also 
Section 4. 



3 General concept of nonlinear self-adjointness 

Motivated by the examples discussed in Sections 1 and 2 as well as other simi- 
lar examples, I suggest here the general concept of nonlinear self-adjointness of 
systems consisting of any number of equations with m dependent variables. This 
concept encapsulates Definition 1.2 of strict self-adjointness and Definition 1.3 
of quasi self-adjointness. The new concept has two different features. They are 
expressed below by two different but equivalent definitions. 

3.1 Two definitions and their equivalence 

Definition 3.1. The system of m differential equations (compare with Eqs. (1-6)) 

= 0, Q; = l,...,m, (3.1) 

with m dependent variables u = (w^, . . . , u^) is said to be nonlinarly self-adjoint 
if the adjoint equations 

F„(x,'u,'y,'U(i),U(i),...,M(,),i;(,)j = — ^-^ = 0, Q; = l,...,m, (3.2) 

are satisfied for all solutions u of the original system (3.1) upon a substitution 

v"' = (p°'{x,u), a = 1, . . . ,m, (3.3) 
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(3.4) 



F*{x,u,ip{x,u), . . . ,u^s),^is)) = Af F^(a;,M, a = l,...,m, 

where are undetermined coefficients, and </7(o.) are derivatives of (3.3), 
^{a) = {Ai ■ ■ ■ Di^{ip"{x,u))}, a = l,...,s. 



(3.5) 



Here v and (p are tlie m-dimensional vectors 



and Eq. (3.4) means tliat not all components (^"(x, u) of vanish simultaneously. 

Remcirk 3.1. If the system (3.1) is over- determined, i.e. fn > m, then the 
adjoint system (3.2) is sub-definite since it contains m <m equations for fn new 
dependent variables v. Vise versa, if m < m, then the system (3.1) is sub-definite 
and the adjoint system (3.2) is over-determined. 

Remark 3.2. The adjoint system (3.2), upon substituting there any solution 
u{x) of Eqs. (3.1), becomes a linear homogeneous system for the new dependent 
variables v°'. The essence of Eqs. (3.5) is that for the self-adjoint system (3.1) 
there exist functions (3.3) that provide a non-trivial (not identically zero) solution 
to the adjoint system (3.2) for all solutions of the original system (3.1). This 
property can be taken as the following alternative definition of the nonlinear 
self- adj ointness . 

Definition 3.2. The system (3.1) is nonlinearly self-adjoint if there exist func- 
tions given by (3.3) that solve the adjoint system (3.2) for all solutions u{x) 
of Eqs. (3.1) and satisfy the condition (3.4). 

Proposition 3.1. The above two definitions are equivalent. 

Proof. Let the system (3.1) be nonlinearly self-adjoint by Definition 3.1. Then, 
according to Remark 3.2, the system (3.1) satisfies the condition of Definition 



Conversely, let the system (3.1) be nonlinearly self-adjoint by Definition 3.2. 

Namely, let the functions f ° given by (3.3) and satisfying the condition (3.4) solve 
the adjoint system (3.2) for all solutions u{x) of Eqs. (3.1). This is possible if 
and only if Eqs. (3.5) hold. Then the system (3.1) is nonlinearly self-adjoint by 
Definition 3.1. 



3.2. 
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Example 3.1. It has been mentioned in Example 1.2 that the KdV equation 

Ut = U^xx + UUx (3.6) 

is strictly self-adjoint. In terms of Definition 3.2 it means that v — u solves the 
adjoint equation 

vt = v^^^ + uvx (3.7) 

for all solutions of the KdV equation (3.6). One can verify that the general 
substitution of the form (3.3), v = ip(t,x,u), satisfying Eq. (3.5) is given by 

V ^ Ai + A2U + As{x + tu), (3.8) 

where Ai, A2, A^ are arbitrary constants. One can also check that v given by Eq. 
(3.8) solves the adjoint equation (3.7) for all solutions u of the KdV equation. 
The solution v = x + tu is an invariant of the Galilean transformation of the KdV 
equation and appears in different approaches (see [10], Section 22.5, and [13]). 
Thus, the KdV equation is nonhnearly self-adjoint with the substitution (3.8). 

Proposition 3.2. Any hnear equation is nonlinearly self-adjoint. 

Proof. This property is the direct consequence of Definition 3.2 because the 
adjoint equation F*[v] = to a linear equation F[u\ — does not involve the 
variable u. 

3.2 Remark on differential substitutions 

One can further extend the concept of self-adjointness by replacing the point-wise 
substitution (3.3) with differential substitutions of the form 

v" = Lp°'{x,u,U(i), . . . ,u^r)), Q; = l,...,m. (3.9) 
Then Eqs. (3.5) will be written, e.g. in the case r = 1, as follows: 

f:{x, u,^,..., u^s), ns)) = Af + A? Dj{F^). (3.10) 
Example 3.2. The reckoning shows that the equation 

Uxy — sinu (3-11) 

is not self-adjoint via a point- wise substitution v — (f{x, y, u), but it is self-adjoint 
in the sense of Definition 3.1 with the following differential substitution: 

V = (p{x, y, Ux, Uy) = Ai[xux - yuy] + A2UX + A^Uy, (3-12) 

where Ai, A2, A3 are arbitrary constants. The adjoint equation to Eq. (3.11) is 

Vxy — V COS U — 0, 

and the self-adjointness condition (3.10) with the function ip given by (3.12) is 
satisfied in the form 

ipxy -ificosu = {Aix + A2)Dx{uxy - sinii) + (A3 - Aiy)Dy{uxy - sinw). (3.13) 



23 



3.3 Nonlinear heat equation 

3.3.1 One-dimensional case 

Let us apply the new viewpoint to the nonlinear heat equation (2.11), Ut — 
{k{u)ux)x, discussed in Section 2.3. We will take it in the expanded form 

ut-k{u)u^^-k\u)ul = 0, k{u)^0. (3.14) 

The adjoint equation (1.18) to Eq. (3.14) is 

Vt + k{u)vxx ^ 0. (3.15) 

We take the substitution (3.3) written together with the necessary derivatives: 

V — (f{t, X, u), 

vt = (fuUt + (ft, = (fuU^ + (fx, (3.16) 
and arrive at the following self-adjointness condition (3.5): 

VuUt + Vt + k{u)[ipuUxx + VuuUl + "^VxuUx + Vxx] 

(3.17) 

= A[ut - k{u)uxx - k'{u)u% 

The comparison of the coefficients of Ut in both sides of Eq. (3.17) yields A = 
Then, comparing the terms with Uxx we see that ipu = 0- Hence Eq. (3.17) reduces 
to 

ipt + k{u)(pxx = (3.18) 

and yields (ft = 0, = 0, whence ip = Cix + C2, where Ci, C2 = const. We 
have demonstrated that Eq. (3.14) is nonlinearly self-adjoint by Definition 3.1 
and that the substitution (3.3) has the form 

v^Cix + C2. (3.19) 

The same result can be easily obtained by using Definition 3.2. We look for 
the solution of the adjoint equation (3.15) in the form v = Lp{t,x). Then Eq. 
(3.15) has the form (3.18). Since it should be satisfied for all solutions u of Eq. 
(3.14), we obtain ip^ = 0, ipxx = 0, and hence Eq. (3.19). 
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3.3.2 Multi-dimensional case 

The similar analysis can be applied to the nonlinear heat equation with several 
variables x — (x^, . . . , x") : 

Ut^V ■ {k{u)Vu) , (3.20) 

or 

Ut-k{u)^u-k\u)\Wu\'^ = Q. (3.21) 
The reckoning shows that the adjoint equation (1.18) to Eq. (3.21) is written 

vt + k{u) Av = 0. (3.22) 

It is easy to verify the nonlinear elf-adjointness by Definition 3.2. Namely, search- 
ing the solution of the adjoint equation (3.22) in the form v — (fi{t,x^, . . . ,x"), 
one obtains 

ifit + k{u) Aip = 0, 

whence 

(ft = 0, Aip = 0. 

We conclude that Eq. (3.21) is self-adjoint and that the substitution (3.3) is given 
by 

= (^(x\...,x"), (3.23) 

where (f{x^, . . . , x") is any solution of the n-dimensional Laplace equation A(p — 
0. 

3.4 Anisotropic nonlinear heat equation 
3.4.1 Two-dimensional case 

Consider the heat diffusion equation 

ut = {f{u)u^)^ + {g{u)uy)y (3.24) 

in an anisotropic two-dimensional medium (see [14], vol. 1, Section 10.8) with 
arbitrary functions f{u) and g{u). The adjoint equation is 

Vt + f{u)v:cx + g{u)vyy = 0. (3.25) 

Using Definition 3.2 we obtain the following equations for nonlinear self-adjointness 
ofEq. (3.24): 

= 0, ip:cx = 0, ipyy = 0. (3.26) 
Integrating Eqs. (3.26) we obtain the following substitution (3.3): 

v^Cixy + C2X + C^y + Ci. (3.27) 
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3.4.2 Three-dimensional Ccise 

The three-dimensional anisotropic nonlinear heat diffusion equation has the fol- 
lowing form (see [14], vol. 1, Section 10.9): 

ut = {f{u)u^)x + {9{u)uy)y + {h{u)u^)^ . (3.28) 

Its adjoint equation is 

vt + f{u)vxx + 9{u)vyy + h{u)v^^ = 0. (3.29) 

Eq. (3.28) is nonlinearly self-adjoint. In this case the substitution (3.27) is 
replaced by 

v^Ci xyz + C2xy + Csxz + C^yz + C5X + Cey + C7 z + Cg. (3.30) 

3.5 Nonlinear wave equations 
3.5.1 One-dimensional case 

Consider the following one-dimensional nonlinear wave equation: 

utt = {k{u)ux)x, k{u) ^ 0, (3.31) 

or in the expanded form 

Utt - k{u)uxx - k'{u)ul = 0. (3.32) 

The adjoint equation (1.18) to Eq. (3.31) is written 

Vtt - k{u)vxx = 0. (3.33) 

Proceeding as in Section 3.3.1 or applying Definition 3.2 to Eqs. (3.32), (3.33) 
by letting v = (f{t,x), wc obtain the following equations that guarantee the 
nonlinear self-adjointncss of Eq. (3.31): 

^tt = 0, ipxx = 0. (3.34) 

Integrating Eqs. (3.34) we obtain the following substitution: 

v^Citx + C2t + Csx + C4. (3.35) 
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3.5.2 Multi-dimensional case 

The multi-dimensional version of Eq. (3.31) with x — (x^, . . . , x^) is written 

Uu^V ■ {k{u)Vu) , (3.36) 

or 

utt - k{u) Au - k'{u)\Vu\'^ = 0. (3.37) 

The adjoint equation is 

vu - k{u) /S.V = 0. (3.38) 

Using Definition 3.2 and searching the solution of the adjoint equation (3.38) in 
the form v — ip{t, x^, . . . , x''), we obtain the equations 

(fitt = 0, = 0. 

Solving them we arrive at the following substitution (3.3): 

V ^ a{x)t + b{x), (3.39) 

where a{x) and b{x) solve the i/-dimensional Laplace equation, 

Aa{x\ ...,x'')^0, Ab{x\ ...,x'')^0. 

Hence Eq. (3.36) is nonlinear ly self- adjoint. 

3.5.3 Nonlinectr vibration of membranes 

Vibrations of a uniform membrane whose tension varies during deformations are 
described by the following Lagrangian: 



L = l 
2 



u^,-k{u)(ul + ul)], k'{u)^0. (3.40) 



The corresponding Euler -Lagrange equation 

dL ^ fdL\ ^ f dL\ ^ f dL 

provides the nonlinear wave equation 

uu = k{u) {u^^ + Uyy) + ^ k\u){ul + ul). (3.41) 

Note that Eq. (3.41) differs from the two-dimensional nonlinear wave equa- 
tion (3.37) by the coefficient 1/2. Let us find out if this difference affects self- 
adjointness. 
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By applying (3.2) to the formal Lagrangian of Eq. (3.41) we obtain: 

1) 

F* = Vtt - k{u) {Vxx + Vyy) - k'{u){UxVx + UyVy+VUxx+VUyy) " " A)" (M ) (M^ + ) . 

We take the substitution (3.3) together with the necessary derivatives (see Eqs. 
(3.16)): 

V ^ ip{t,x,y,u), vt^ (puUt + (ft, 

Vx = ^uUx + VX, Vy = (PuUy + Ify, 

Vxx = VuUxx + ^uuUl + 'iipxuUx + ^xx, (3.42) 
^yy — ^u^yy ~l~ ^utL^y ~l~ ^^yu^y ~l~ ^yy^ 
Vtt = ^uUtt + VuuUl + '^VtuUt + ^tU 

and substitute the expressions (3.42) in the self-adjointness condition (3.5): 

F\=^ = >\utt - k{u) {uxx + Uyy) - ^ k'{u){ul + ul)]. 

Comparing the coefficients of Utt we obtain X = Lp^- Then wc compare the co- 
efficients of Uxx and obtain ipk'{u) = 0. This equation yields Lp = because 
k'{u) 7^ 0. Thus, the condition (3.4) is not satisfied for the point-wise substitu- 
tion (3.3). Further investigation of Eq. (3.41) for the nonlinear self-adjointness 
requires differential substitutions. 

3.6 Anisotropic nonlinear wave equation 
3.6.1 Two-dimensional case 

The two-dimensional anisotropic nonlinear wave equation is (see [14], vol. 1, 
Section 12.6) 

Utt = {f{u)ux)x + {g{u)uy)y . (3.43) 
Its adjoint equation has the form 

Vtt - f{u)vxx - g{u)vyy = 0. (3.44) 

Proceeding as in Section 3.4 we obtain the following equations that guarantee the 
self-adjointness of Eq. (3.43): 

Vtt = 0, (fxx = 0, (fyy = 0. (3.45) 

Integrating Eqs. (3.45) we obtain the following substitution (3.3): 

v^Ci txy + C2tx + Caty + C^xy + C5t + Cqx + C7y + Cg. (3.46) 
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Remcirk 3.3. I provide here detailed calculations in integrating Eqs. (3.45). 
The general solution to the linear second-order equation iptt = is given by 

ip^A{x,y)t + B{x,y) (3.47) 

with arbitrary functions A{x,y) and B{x,y). Substituting this expression for (f 
in the second and third equations (3.45) and splitting with respect to t we obtain 
the following equations for A{x,y) and B{x,y) : 

A _ fi A _ fi 

^xx — ^yy — 
Substituting the general solution 

A = ai{y)x + a2{y) 

of the equation A^^ = in the equation Ayy = and splitting with respect to x, 
we obtain a" — 0, 03 = 0, whence 

ai = cny + C12, 02 = C2iy + C22, 

where cn, . . . , C22 are arbitrary constants. Substituting these in the above expres- 
sion for A we obtain 

A = ciixy + C12X + C2iy + C22. 
Proceeding likewise with the equations for B{x,y), we have 

B = duxy + di2X + d2iy -\- ^22 

with arbitrary constant coefficients cin, . . . , ^22- Finally, we substitute the result- 
ing A and B in the expression (3.47) for (/? and, changing the notation, arrive at 
(3.46). 

3.6.2 Three-dimensional case 

The three-dimensional anisotropic nonlinear wave equation 

'f^tt = {f{u)ux)x + {g{u)uy)y + {h{u)uz)z (3.48) 
has the following adjoint equation 

vtt - f{u)v^^ - g{u)vyy - h{u)Vzz = 0. (3.49) 
In this case Eqs. (3.45) are replaced by 

and yield the following substitution (3.3): 

V — Ci txyz + C2 txy + C3 txz + tyz + C^tx + Cq ty + C7 tz 

(3.50) 

+ C8xy + Cg xz + Clo yz + Cn t + Cu x + C13 y + Cu z + C15. 



29 



3.7 Nonlinear self-adjointness and multipliers 



The approach of this section is not used for constructing conservation laws. But 
it may be useful for other applications of the nonlinear self-adjointness. 



Theorem 3.1. The differential equation (1.17), 

F{x,u,u^i),...,u^s)) = 0, 



(3.51) 



is nonlinearly self-adjoint (Definition 3.1) if and only if it becomes strictly self- 
adjoint (Definition 1.2) upon rewriting in the equivalent form 



IJ,{X, u)F{x, U, 14(1), ■ ■ ■ , U^s)) = 0, iJ,{x, u) ^ 0, 

with an appropriate multiplier ij,{x,u). 



(3.52) 



Proof. We will write the condition (3.5) for nonlinear self-adjointness of Eq. 
(3.51) in the form 



S{vF) 



Su 



v=<p{x,u) 



(3.53) 



Furthermore, invoking that the equations (3.52) and (3.51) are equivalent, we 
will write the condition (1.28) for strict self-adjointness of Eq. (3.52) in the form 



6{wfiF) 



6u 



= X{x,u)F{x,U,U(^1), . . .,U(^s))- 



(3.54) 



Since -u; is a dependent variable and fi = fi{x,u) is a certain function of x,u, the 
variational derivative in the left-hand side of (3.54) can be written as follows: 



5{wiiF) dfj, 



dF 



5u 



du 

w^F + 
ou 



du 

S(vF) 
Su 



dF\ 



w— F + iiw— Di I iiw— — -|- DiDj I ixw 



duij 



OF 



du 



where v is the new dependent variable defined by 

V — ii{x, u)w. 



(3.55) 



is the new dependent variable instead of w. Now the left side of Eq. (3.54) is 
written 



6{wi2F) 



5u 



du 6u 



(3.56) 



v=uiJ.{x,u) 
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Let us assume that Eq. (3.51) is nonlinearly self-adjoint. Then Eq. (3.53) 
holds with a certain given function (p{x,u). Therefore, we take the multiplier 

Mx,«) = ^ (3.57) 
u 

and reduce Eq. (3.56) to the following form: 



5u 

This proves that Eq. (3.54) holds with 



ou u 



OU u 

Hence, Eq. (3.52) with the multiplier ^ given by (3.57) is strictly self-adjoint. 

Let us assume now that Eq. (3.52) with a certain multiplier //(x, u) is strictly 
self-adjoint. Then Eq. (3.54) holds. Therefore, if we take the function (/? defined 
by (see (3.57)) 

ip{x,u) = uiJ,{x,u), (3.58) 

Eq. (3.56) yields: 

S{vF) 



v=ip{x,u) 

It follows that Eq. (3.53) holds with 



ou 



A = A - u—. 

ou 

We conclude that Eq. (3.51) is nonhnearly self-adjoint, thus completing the proof. 

Example 3.3. The multiplier (2.22) used in Example 2.3 and the function ip — 
x'^ that provides a solution of the adjoint equation (2.21) to the Kompaneets 
equation are related by Eq. (3.58). 

Example 3.4. Let us consider the one-dimensional nonhnear wave equation 
(3.32), 

Utt - k{u)u^^ - k'{u)ul = 0. 

If we substitute in (3.57) the function (p given by the right-hand side of (3.35) 
we will obtain the multiplier that maps Eq. (3.32) into the strictly self-adjoint 
equivalent form. For example, taking (3.35) withCi = C3 = C4 = 0, C2 = 1 we 
obtain the multiplier 

t 

-• 

u 
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The corresponding equivalent equation to Eq. (3.32) has the formal Lagrangian 



tv 

C = —[utt - k{u)uxx - k'{u)u 
u 



We have 
5C 



5u ^'1 



ftv 



tv 



— k\u) k [u)Uxx + — k{u)Ux 



+ 2Dx {—k'{u)uA - — k"{u)ul + % k\u)ul. 
\ u J u 



Letting here v = u we see that the strict self-adjointness condition is satisfied in 
the following form: 



5£ 

5u 



--[utt - k{u)uxx - k'{u)u 
u 



4 Generalized Kompaneets equation 

4.1 Introduction 

The equation 



aI dn 2 
X I — + n + n 

ox 



(4.1) 



dn Id 
dt x"^ dx 

known as the Kompaneets equation or the photon diffusion equation, was derived 
independently by A.S. Kompaneets^ [15] and R. Weymann [16]. They take as 
a starting point the kinetic equations for the distribution function of a photon 
gas^ and arrive, at certain idealized conditions, at Equation (4.1). This equation 
provides a mathematical model for describing the time development of the energy 
spectrum of a low energy homogeneous photon gas interacting with a rarefied 
electron gas via the Compton scattering. Here n is the density of the photon gas 
(photon number density) , t is time and x is connected with the photon frequency 
V by the formula 

hi' , , 

where h is Planck's constant and kT^. is the electron temperature with the standard 
notation k for Boltzmann's constant. According to this notation, hv has the 



^He mentions in his paper that the work has been done in 1950 and pubhshed in Report N. 
336 of the Institute of Chemical Physics of the USSR Acad. Sci. 

^Weymann uses Dreicer's kinetic equation [17] for a photon gas interacting with a plasma 
which is slightly different from the equation used by Kompaneets. 
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meaning of the photon energy. The nonrelativistic approximation is used, i.e. 
it is assumed that the electron temperatures satisfy the condition kT^. ^ mc^, 
where m is the electron mass and c is the light velocity. The term low energy 
photon gas means that hu <^ mc^. 

The question arises if the idealized conditions assumed in deriving Eq. (4.1) 
may be satisfied in the real world. For discussions of theoretical and observa- 
tional evidences for such possibility in astrophysical environments, for example 
in intergalactic gas, see e.g. [18], [19] and the references therein. See also the 
recent publication 



4.2 Discussion of self-adjoint ness of the Kompaneets equa- 
tion 

For unifying the notation, the dependent variable n in Eq. (4.1) will be denoted 
by u and Eq. (4.1) will be written further in the form 

ut = ^D^lx^u^ + u + u^)]. (4.3) 
Writing it in the expanded form 

Ut = x'^u^^c + {x'^ + 4:X + 2x^u)ux + 4:x{u + u^), (4.4) 
we have the following formal Lagrangian for Eq. (4.3): 

C — v[—Ut + x^Uxx + {x^ + 4a; + 2x^u)ux + ^x{u + u^)]- 
Working out the variational derivative of this formal Lagrangian, 

— = Dt{v) + Dl{x^v) - Dx[{x^ + 4x + 2x^u)v\ + 2x^vUx + 4x(l + 2u)v, 

(J Lv 

we obtain the adjoint equation to Eq. (4.3): 
SjC 

— = vt + x^Vxx - + 2u)vx + 1{x + 2xu - \)v = 0. (4.5) 
\iv — ip{u)^ then 

vt = ^'{u)ut, Vx = (p'{u)Ux, Vxx = v'{u)uxx + (p"{u)ul. 
It follows that the quasi self-adjointness condition (1.36), 

= \[—Ut + x'^Uxx + {x^ + 4x + 2x^u)ux + 'ix{u + w^)]. 



5_C 

Su 



v=ip(u) 
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is not satisfied. 

Let us check if this condition is satisfied in the more general form (3.5): 

— ^ X[-Ut + x^u^^ + {x^ + Ax + 2x^u)u:, + Ax{u + u^)]. (4.6) 

v=ip{t,x,u) 

In this case 

Vt = Dt[ip{t, X, u)] = (puUt + ipt, 

Vx = D^[(p{t, X, U)] = (fiuUx + <Px: (4.7) 
'^xx — DxiVx) — 'fiu'^xx ~l~ '^uu^x '^'^xu^x ~l~ '^xx- 

Inserting (4.7) in the expression for the variational derivative given by (4.5) and 
singling out in Eq. (4.6) the terms containing ut and Uxx, we obtain the following 
equation: 

^u{ut + X^Uxx] = M-Ut + x'^Uxx]- 

Since this equation should be satisfied identically in Ut and Uxx, it yields A = 
ipu = 0. Hence ip = (p{t,x) and Eq. (4.6) becomes: 

ipt + x^ifxx - a^^(l + 2m)v5:« + 2(x + 2xu - !)(/? = 0. (4.8) 

This equation should be satisfied identically in x and u. Therefore we nullify 
the coefficient for u and obtain 

- 0, 

whence 

V9(t, X) = c{t)x'^. 

Substitution in Eq. (4.8) yields c'{t) = 0. Hence, v = (p{t, x) = Cx'^ with arbitrary 
constant C. Since A = in (4.6) and the adjoint equation (4.5) is linear and 
homogeneous in v, one can let C = 1. Thus, we have demonstrated the following 
statement. 

Proposition 4.1. The adjoint equation (4.5) has the solution 

V = x^ (4.9) 

for any solution u of Equation (4.3). In another words, the Kompaneets equation 
(4.3) is nonlinearly self-adjoint with the substitution (3.3) given by (4.9). 

Remark 4.1. The substitution (4.9) does not depend on u. The question arises 
on existence of a substitution v = ip{t, x, u) involving u if we rewrite Eq. (4.3) in 
an equivalent form 

a{t, X, u) [-Ut + x^Uxx + {x^ + 4a; + 2x^u)ux + Ax{u + u^)] = (4.4') 

with an appropriate multiplier a 7^ 0. This question is investigated in next section 
for a more general model. 
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4.3 The generalized model 

In the original derivation of Eq. (4.1) the following more general equation appears 
accidentally (see [15], Eqs. (9), (10) and their discussion): 



dn Id 



dt g{x) dx 



/(a:)(g + /(n; 



(4.10) 



with undetermined functions f{u) and g{x). Then, using a physical reasoning, 
Kompaneets takes f{u) = n(l + n) and g{x) = . This choice restricts the 
symmetry properties of the model significantly. Namely, Equation (4.1) has only 
the time-translational symmetry with the generator 

The symmetry (4.11) provides only one invariant solution, namely the stationary 
solution n — n{x) defined by the Riccati equation 

dn C 
— + n +n = —r ■ 
ax X* 

The generalized model (4.10) can be used for extensions of symmetry properties 
via the methods of preliminary group classification [21, 22]. In this way, exact 
solutions known for particular approximations to the Kompaneets equation can 
be obtained. This may also lead to new approximations of the solutions by taking 
into account various inevitable perturbations of the idealized situation assumed 
in the Kompaneets model (4.1). 

So, we will take with minor changes in notation the generalized model (4.10): 

'^t = j^ D,{h^{x)[u, + f{u)]}, h\x) ^ 0. (4.12) 

It is written in the expanded form as follows: 

ut = h{x) («,, + f{u)u,) + 2h'{x) {u, + f{u)) . (4.13) 
We will write Eq. (4.13) in the equivalent form similar to (4.4'): 

a(t, x/a)[-ut + h{x) {u^^ + f{u)u^) + 2h\x) {u^ + /(«))] = 0, (4.14) 
where a 7^ 0. This provides the following formal Lagrangian: 

C = v a{t, x,u)[-ut + h{x) [u^^ + f'{u)u^) + 2h'{x) [u^ + /(«))] , (4.15) 
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= X[-Ut + (u,, + fu,)h + 2(u, + f)h'] (4.17) 

V=ip{t,X,U) 



where v is a new dependent variable. For this Lagrangian, we have 

— = Dt(va) + Dl[h(x)va] - DJh(x)f(u)va + 2h'(x)va] 
du 

+ h{x)f"{u) vau:, + 2h\x)f{u)va 

+ vau[- ut + h{x) {u:cx + f{u)u^) + 2h'{x) {u^ + /(«))] . 
The reckoning shows that 

r p 

— = Dtiva) + hDliva) - hfD^{va) + {h'f - h")va 
du 

+ vau[-ut + (u^^ + f'u^)h + 2{u^ + f)h']. (4.16) 
Now we write the condition for the self-adjointness of Eq. (4.13) in the form 

6C 
du 

with an undetermined coefficient A. Substituting (4.16) in (4.17) we have: 

Dtiipa) + hDli^a) - hfD^{^a) + {h'f - h")ipa 
+ ipau[-ut+ {u^^ + fu^)h + 2{u^ + f)h'] (4.18) 
^\[-Ut + {u^^ + fu^)h + 2{u^ + f)h']. 
Here ip = ip{t,x,u), a = a(t,x,u) and consequently (see (4.7)) 

A(v^tt) = {'^a)uUt + (v?a)t, 

D^{ipa) = {ipa)uU^ + {(pa)x, (4.19) 
Dl{cpa) = {(pa)uU^^ + {(pa)uuul + 2{(pa)^uUx + {.^ot)^^. 

We substitute (4.16) in Eq. (4.18), equate the coefficients for ut in both sides of 
the resulting equation and obtain {(pa)u — ^ocu — —A. Hence, 

A = 

Using this expression for A and equating the coefficients for hu^x in in both sides 
of Eq. (4.18) we get (^^a)^ + = —a^u- It follows that {(pa.)u — and hence 

aip — k{t, x). 

Now Eq. (4.18) becomes: 

kt + h{x)kxx - h'\x)k + f\u)[h\x)k - h{x)kx] = 0. 
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If f"{u) 7^ 0, the above equation splits into two equations: 

h'{x)k — h{x)kx = 0, kt + h{x)kxx — h"{x)k. 

The first of these equations yields k{t, x) = c{t)h{x), and then the second equation 
shows that c'{t) = 0. Hence, k = C h{x) with C = const. Letting C = 1, we have: 



a(f — h{x). 
Eq. (4.20) can be satisfied by taking, e.g. 

h{x) 



a 



u 



if = u. 



(4.20) 



(4.21) 



Thus, we have proved the following statement. 
Proposition 4.2. Eq. (4.12) written in the equivalent form 

^Ut^-D,{h\x)[u, + f{u)]] 

Ui Hi 



(4.22) 



is strictly self-adjoint. In another words, the adjoint equation to Eq. (4.22) 
coincides with (4.22) upon the substitution 



V = u. 



(4.23) 



In particular, let us verify by direct calculations that the original equation 
(4.3) becomes strictly self-adjoint if we rewrite it in the equivalent form 

(4.24) 



u u ^ 



Eq. (4.24) reads 



x^ 

Ut-\ Uxx + 

u u 



(x^ + Ax^)- + 2x* Ux + 4x^(1 + -u) = 
u J 



(4.25) 



and has the formal Lagrangian 



2V 4U 

— — X — Ui ~\~ X — "^xx ~t~ 

u u 



(x^ + Ax^) - + 2x% Ux + Ax^{v + uv). 

XL 



Accordingly, the adjoint equation to Eq. (4.25) is written 



ix^ + Ax'')- + 2x 
u 



+ X^^Ut -X^^Uxx - {x^ + Ax^)^Ux + Ax^v = 0. 



Letting here v = u one has v/u = 1 and after simple calculations arrives at Eq. 
(4.25). 
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5 Quasi self-adjoint reaction-diffusion models 



Let us consider the one-dimensional reaction-diffusion model described by the 
following system (see, e.g. [23]): 



(5.1) 



(5.2) 



(5.3) 



- = nu,.)+A—^+-(^^M-y 

0^-9M + B— + -i^^^iu,v)-j. 

It is convenient to write Eqs. (5.1) in the form 

Dtiu) - ADliu) + [<f>{u, v)D,{v)] + f{u, v), 

Dt{v) = BDliv) + D, [iPiu, u)DM] + 9{u. 
The total differentiations have the form 

^ d d d d d d d 

Dt^ — + Ut^ + Vt^ + Utt^ h Utx^ h Vtt^ + Vtx^ \ , 

Ot OU OV OUt OUx OVt OVx 

d d d d d d d 

Dx — t; \- Ux-^ 1" VxT^ — h Utx-^ 1" tJ-xx-^ 1" 't'txT; 1" '^xx-^ 1" ■ ■ 

ox OU OV OUt OUx OVt (JVx 

and Eqs. (5.2) are written 

Ut = AUxx + (f>Vxx + [(f>uUx + (pvVx] Vx + f, 
Vt = BVxx + i^Uxx + [i^uUx + i^vVx] Ux + 9- 

The formal Lagrangian for the system (5.4) is 

C = z{AUxx -Ut + (j)Vxx + (t>uUxVx + (t>vVl + /) 
+ w{BVxx -Vt + l/jUxx + l/^uUl + ll^vUxVx + 5'), 

where z and w are new dependent variables. Eqs. (1.8) are written: 

Su ^""^duxx) ^^^dut) ^''^dux)^du' 

^ 5v ^ \dvxxj * \dvtj ^ \dvxj dv 
Substituting here the expression (5.5) for C we obtain after simple calculations 
the following adjoint equations (3.2) to the system (5.4): 

AZxx + Zt + l/jyVxWx - (puVxZx + ll^Wxx + Zfu + WQ^ = 0, (5.6) 
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(5.4) 



(5.5) 



BWxx + Wt + (puUxZx - IpvUxWx + (t>Zxx + zfy + WQy = 0. (5.7) 

Let us investigate the system (5.4) for quasi self-adjointness (Definition 1.3). 
We write the left-hand sides of Eqs. (5.6) and (5.7) as hnear combinations of the 
left-hand sides of Eqs. (5.4): 

AZxx + Zt + IpvVxWx - (puVxZx + IpWxx + zfu + WQu (5.8) 

= {AUxx -Ut + (pVxx + (puUxVx + (pvVl + f)P 

+ {Bv^^ -vt + ijjuxx + i>uul + i>vUxVx + g)Q, 

BWx:x + Wt + (I) uUxZx - l/jyUxWx + (pZxx + Zfy + WQy (5.9) 
= {AUxx -Ut + (t)Vxx + (t>uUxVx + (t>vVl + f)M 
+ {BVxx -Vt + IpUxx + i^uUl + i^vUxVx + g)N, 

where P, Q, M and N are undetermined coefficients. We write the substitution 
(1.42) in the form 

z = Z{u, v), W{u, v) (5.10) 

and insert in the left-hand sides of Eqs. (5.8)-(5.9) these expressions for z,w 
together with their derivatives 

Zt = ZuUt + ZyVti = ZuUx + ZyVxi 

ZxX Zy^UxX ~l~ ZyVxX ~l~ Z^yU^ -\- 2,ZyyUxVx ~\~ ZyyV^, 

Wt = WuUt + WyVt, Wx = WuUx + WyVx, 

WxX = WuUxx + WyVxX + WuuUl + 2WuyUxVx + WyyVj. 

Equating the coefficients for Uf and Vt in both sides of Eqs. (5. 8)- (5. 9) we obtain 

P — ~Zu, Q — —Zy, /rii\ 

5.11) 

N^-Wy, M^-Wy. ^ ' 

Now we calculate the coefficients for Uxx and Vxxi take into account Eqs. (5.11) 
and arrive at the following equations: 

2AZy + ^l;Zy + i;Wy = 0, 

{A + B)Zy + ^Zy + i^Wy = 0, 

(5.12) 

2BWy + (t)Zy + i:Wy = 0, 

{A + B)Wu + 4>Z^ + ^Wy = Q. 
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Eqs. (5.12) provide a linear homogeneous algebraic equations for the quanti- 
ties 

Zu, Z,, Wu, Wu 

with the matrix 

(2A V ^ \ 

A + B tl) 

(j) 2B ' 

A + B J 

This matrix has the inverse because its determinant is equal to 

A{A + Bf{(f)ijj-AB) 

and does not vanish in the case of arbitrary A,B,4> and ip. Hence, Eqs. (5.12) 
yield: 

Z^^Z.^Wu^W^^ 0. (5.13) 

It follows that Z{u, v) = Ci, W{u, v) = C2. Thus, the substitution (1.42) has the 
form 

z^Ci, w^C2 (5.14) 
with arbitrary constants Ci, C2. Then Eqs. (5.8)-(5.9) become 

(Ci/ + C2g)u = 0, (Ci/ + C2g), = 

and yield 

where / = Ci/, g = C25', and C = const. Since / and g, along with / and g, are 
arbitrary functions, we can omit the "tilde" and write 

f + g^C. (5.15) 

Eq. (5.15) provides the necessary and sufficient condition for the quasi self- 
adjointness of the system (5.1). Thus, we have proved the following statement. 

Theorem 5.1. The system (5.1) is quasi self-adjoint if and only if it has the 
form 

^ = /(-,-) + ^^ + ^(0(-,-)^j, 

(5.16) 

where (f){u,v), iIj{u,v), f{u,v) are arbitrary functions and A,B,C are arbitrary 
constants. The substitution (1.42) is given by (5.14). 
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Remcirk 5.1. If we replace (5.10) by the general substitution (3.3), i.e. take 



z = Z{t,x,u,v), w — W{t,x,u,v), 
then Eqs. (5.14) will be replaced by 

z^Z{t,x), w^W{t,x), 
with functions Z{t,x), W{t,x) satisfying the following equations: 



(5.17) 



(5.18) 



{AW - <PuZ)^ = 0, 

AZ,, + Zt + Wxx + {fZ + gW)u = 0, 



(5.19) 



(5.20) 



BW,, + Wt + ct>Z,:, + {fZ + gW) 



V 



0. 



6 A model of an irrigation system 



Let us consider the second-order nonlinear partial differential equation 



(6.1) 



It serves as a mathematical model for investigating certain irrigation systems 
(see [14], vol. 2, Section 9.8 and the references therein). The dependent variable 
ip denotes the soil moisture pressure head, C{ip) is the specific water capacity, 
K{iIj) is the unsaturated hydraulic conductivity, S^ip) is a source term. The 
independent variables are the time t, the horizontal axis x and the vertical axis 
z which is taken to be positive downward. 

The adjoint equation (3.2) to Eq. (6.1) has the form 



It follows from (6.2) that Eq. (6.1) is not nonlincarly self-adjoint if C (ip) , K {ip) 
and S^ip) are arbitrary functions. Indeed, using Definition 3.2 of the nonlinear 
self-adjointness and nullifying in (6.2) the term with S'{ijj) we obtain v = 0. 
Hence, the condition (3.4) of the nonlinear self-adjointness is not satisfied. 

However, Eq. (6.1) can be nonlinearly self-adjoint if there are certain relations 
between the functions C(V'), K{i/j) and S^ip). For example, let us suppose that the 
specific water capacity C^ip) and the hydraulic conductivity K^ip) are arbitrary, 
but the source term S{ip) is related with C('V') by the equation 



(6.2) 



S\'4>) — aC{ip), a — const. 



(6.3) 



Then Eq. (6.2) becomes 

cmvt 



av] + K{ilj) [v^^ + V,,] + K'{ilj)v, = 
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and yields: 

Vz = 0, = 0, vt-av^ 0. (6.4) 
We solve the first two equations (6.4) and obtain 

V = p{t)x + q{t). 

We substitute this in the third equation (6.4), 

\p\t) - ap{t)]x + q'it) - aq{t) = 0, 

split it with respect to x and obtain: 

p{t) - apit) = 0, q{t) - aq{t) = 0, 

whence 

p(t) = q{t) = /e"*, b, I = const. 

Thus, Eq. (6.1) satisfying the condition Eq. (6.3) is nonlinearly self-adjoint, and 
the substitution (3.3) has the form 

v^(bx + Oe"*. (6.5) 

One can obtain various nonlinearly self-adjoint Equations (6.1) by considering 
other relations between C (ijj) , K (ijj) and <S'('0) different from (6.3). 
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Part 2 

Construction of conservation laws 
using symmetries 

7 Discussion of the operator identity 

7.1 Operator identity and alternative proof of Noether's 
theorem 

Let us discuss some consequences of the operator identity^ 

x + A(r) = w^"^ + ANV (7.1) 



Here 



W = v'-eu^, a = l,...,m, (7.3) 
- + 5^(-irA,---A.^-^, « = l,...,m, (7.4) 



SW^ dw^ ' ' dul , 



and 



where the Euler-Lagrange operators with respect to derivatives of are obtained 
from (7.4) by replacing by the corresponding derivatives, e.g. 

S^ = e^ + YS-^fD„..D, ■ (7.6) 

* * S=l '■Jl Js 

The coefficients in (7.2) arc arbitrary differential functions (see Section 

1.1) and the other coefficients are determined by the prolongation formulae 

cr = D,(w^) + e<^, c^,^ = A, A.(w^") + eu%,,„ .... (7.7) 



■^Recently I learned that the identity (7.1) was proved in [24]. Namely, Eq. (7.1) is the same 
(except for notation) as Eq. (19) from [24]. The operator identity (7.1) was rediscovered in [25] 
and used for simplifying the proof of Noether's theorem. Accordingly, Eq. (7.1) was called in 
[25] the Noether identity. See also [4], Section 8.4. 
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The derivation of Eq. (7.1) is essentially based on Eqs. (7.7). 

Recall that Noether's theorem, associating conservation laws with symme- 
tries of differential equations obtained from variational principles, was originally 
proved by calculus of variations. The alternative proof of this theorem given in 
[25] (see also [10, 4]) is based on the identity (7.1) and is simple. Namely, let us 
consider the Euler- Lagrange equations 



5C 



0, a 



, m. 



(7.8) 



If we assume that the operator (7.2) is admitted by Eqs. (7.8) and that the 
variational integral 

is invariant under the transformations of the group with the generator X then 
the following equation holds: 

X(£) + A(e)>C = 0. (7.9) 
Therefore, if we act on C by both sides of the identity (7.1), 

X{C) + D,{C)C = + D.mC)] , 

and take into account Eqs. (7.8), (7.9), we see that the vector with the compo- 
nents 

C^ = N'(£), i = l,...,n, (7.10) 
satisfies the conservation equation 



(7.8) 



0. 



(7.11) 



For practical applications, when we deal with law order Lagrangians C, it is 
convenient to restrict the operator (7.5) on the derivatives involved in C and 
write the expressions (7.10) in the expanded form 



dC 



dC 



du 



+ D.Dk 



dC 



(7.12) 



dC 
dufj 



dC 
9Kjk 



+ ... 



dC 



Thus, Noether's theorem can be formulated as follows. 
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Theorem 7.1. If the operator (7.2) is admitted by Eqs. (7.8) and satisfies the 
condition (7.9) of the invariance of the variational integral, then the vector (7.12) 
constructed by Eqs. (7.12) satisfies the conservation law (7.11). 

Remark 7.1. The identity (7.1) is valid also in the case when the coefficients 
?7" of the operator X involve not only the local variables x, m, «,(!), M(2), .. . 
but also nonlocal variables (see Section 11.5). Accordingly, the formula (7.12) 
associates conserved vectors with nonlocal symmetries as well. 

Remark 7.2. If the invariance condition (7.9) is replaced by the divergence con- 
dition 

X{jr) + D,{C)C^D,{B'), 

then the identity (7.1) leads to the conservation law (7.11) where the conserved 
vector (7.10) is replaced with 

C' = N\j:) - B\ i = l,...,n. (7.13) 

Remark 7.3. If we write the operator (7.2) in the equivalent form 

then the prolongation formulae (7.7) become simpler: 

Cr = D,{W"), Q, = AiA.(w-"),.... (7.15) 

7.2 Test for total derivative and for for divergence 

I recall here the well-known necessary and sufficient condition for a differential 
function to be divergence, or total derivative in the case of one independent 
variable. 

One can easily derive from the definition (1.1) of the total differentiation Di 
the following lemmas (see also [4], Section 8.4.1). 

Lemma 7.1. The following infinite series of equations hold: 

^ D. = D. ^ 



a 



du"' du 



D.DkJ^ A = A A + AA-Dfc ^ 
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Lemma 7.2. The following operator identity holds for every i and a : 

S 



A = 0. 



5m" 

Proof. Using Lemma 7.1 and manipulatinng with summation indices we obtain: 

S f d d d d \ 
-^Di = D.- h A-DfeT^ D.DkDi- h • • • A 

d d d d d 
= -T^Di - Di DiD .—- + DiDk ^ h DiD.Dk^ — 

-DiDkDi^----^0. 

Proposition 7.1. A differential function f{x, u, U(i), . . . , Uf^s)) & A\s divergence, 

f = Di{h'), h\x,u,...,u^s-i))^A, (7.16) 
if and only if the following equations hold identically in x, u, . . . : 

Sf 



= 0, a^l,...,m. (7.17) 



The statement that (7.16) implies (7.17) follows immediately from Lemma 
7.2. For the proof of the inverse statement that (7.17) implies (7.16), see [26], 
Chapter 4, § 3.5, and [24]. See also [4], Section 8.4.1. 

We will use Proposition 7.1 also in the particular case of one independent 
variable x and one dependent variable u — y. Then it is formulated as follows. 

Proposition 7.2. A differential function f{x,y,y', . . . ,y''^^) G ^ is the total 
derivative, 

/ = A(y), gix,y,y',...,y^'-'^)eA, (7.18) 
if and only if the following equation holds identically in x,y,y', .. . : 

= 0. (7.19) 

dy 

Here Sf/6y is the Euler- Lagrange operator (7.6): 

Sy~dy ''dy'^ ^dy" ^ dy'" ^ ' ^ ^ ^ 
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7.3 Adjoint equation to linear ODE 

Let us consider an arbitrary sth-order linear ordinary differential operator 

L[y] = aoy^'^ + a^y^'-'^ + ■■■ + a,_2|/" + a,_i2/' + a,y, (7.21) 

where Oj = ai{x). We know from Section 1.5 that the adjoint operator to (7.21) 
can be calculated by using Eq. (1.8). I give here the independent proof based on 
the operator identity (7.1). 

Proposition 7.3. The adjoint operator to (7.21) can be calculated by the for- 
mula 

L*[.] = ^^- (7.22) 
dy 

Proof. Let 

be the operator (7.14) with one independent variable x and one dependent vari- 
able u = y, where the prolongation formulae (7.15) are written using the notation 

w'^D,{w), w"^Dl{w),.... (7.24) 

In this notation the operator (7.5) is written 

Having in mind its application to the differential function L[y] given by (7.21) we 
consider the following restricted form of N : 

N = w/- + «;'/- + ■■■ + . (7.25) 

Sy' Sy" Sy('^ ^ ' 

The identity (7.1) has the form 

r 

X^w—^D^n. (7.26) 
by 

We act by both sides of this identity on zL\y\ where 2; is a new dependent variable: 

X{zL\y\) = w^-^^^ + DM^L[y]) . (7.27) 
dy 

Since the operator (7.23) does not act on the variables x and z, we have 

XizL[y]) = zXiL[y]). (7.28) 
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Furthermore we note that 

X{L[y]) ^ L[w]. 
Inserting (7.28) and (7.29) in Eq. (7.27) we obtain 

5{zL[y]) 



zL[w\ — w- 



5y 



where is a quadratic form ^ — ^[w,z] defined by 

^ = N{zL[y]). 



(7.29) 



(7.30) 



(7.31) 



After replacing w with y Eq. (7.30) coincides with Eq. (1.3) for the adjoint 
operator, 

zL[y]-yL*[z] = D,{ij), (7.32) 
where L*[z] is given by the formula (7.22) and ip — ip[y, z] is defined by 



^P[y,z] = ^w,z]l^^ = HizL[y])l^^. 



(7.33) 



Remark 7.4. Let us find the explicit formula for ijj in Eq. (7.32) We write the 
operator N given by Eq. (7.25) in the expanded form 



N = w 

+ w' 
+ w 



dy' ^dy"^ dy^^) 
' ' -D.^ + --- + (-D.)- ' 



{s-2) 



dy" ^ dy 

d 



d 



dyi') 



+ w 



{s-l) 



dy(^) 

5 

' 



+ 



_dy('-^^ 

act on zL[y] written in the form 

zL[y] = Usyz + ag-iy'z + as-2y"z H h aiy'^^~'^h + aoy'^^h, 

and obtain ^. We replace w with y in — "i/lw, z] and ip = ip[y,z] : 



'^-\aozf-'^] 



ip[y, z]^y [us-i z - {as-2 ^)' H h (-1^ 

+ y' [as-2 z - (a,_3 z)' + ■ ■ ■ + (-l)^-'(ao z)^'-'^] + 
+ yM [a,z-{aozy]+y^'-'^ a^z. 



(7.34) 



The expression (7.34) is obtained in the classical literature using integration 
by parts (see, e.g. [27], Chapter 5, §4, Eq. (31')). 
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7.4 Conservation laws and integrating factors for linear 
ODEs 

Consider an sth-order homogeneous linear ordinary differential equation 

L[y] = 0, (7.35) 
where L[y\ is the operator defined by Eq. (7.21). If L[y\ is a total derivative, 



Eq. (7.35) is written as a conservation law 

(^(x,?/,y',...,y(^-^)))=0, 
whence upon integration one obtains a linear equation of order s — 1 : 

V'(x,y,y',...,y(^-^)) = Ci. 
We can also reduce the order of the non-homogeneous equation 

m = f{x) 

by rewriting it in the the conservation form 

V'(x, y,y',..., y^'-^'^) - j f{x)dx 



(7.36) 



0. 



(7.37) 



(7.38) 



(7.39) 



Integrating it once we obtain the non-homogeneous linear equation of order s — 1 
^l:{x,y,y\...J'''^) = C,+ f f{x)dx. 



Example 7.1. Consider the second-order equation 

y" + y' sinx + y cos x — 0. 

We have 

y" + y' sinx + ycosx = D^iy' + ysmx). 

Therefore the second-order equation in question reduces to the first-order equa- 
tion 

y' + ysvax = Ci. 
Integrating the latter equation we obtain the general solution 



y 



C2 + C1J e 
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' dx 



to our second-order equation. Dealing likewise with the non-homogeneous equa- 
tion 

7/" + y' sin X + y cos x = 2x 
we obtain its general solution 



y 



C2 + J {Ci + x^) e 



dx 



If L[y] in Eq. (7.35) is not a total derivative, one can find an appropriate 
factor (f){x) 7^ 0, called an integrating factor, such that (f){x)L[y] becomes a total 
derivative: 

mm = {^{x, y,y',..., y(^-i))) . (7.40) 

A connection between integrating factors and the adjoint equations for linear 
equations is well known in the classical literature (see, e.g. [27], Chapter 5, §4). 
Proposition 7.2 gives a simple way to establish this connection and prove the 
following statement. 

Proposition 7.4. A function (f)(x) is an integrating factor for Eq. (7.35) if and 
only if 

z^cj>{x), 0(x)^O, (7.41) 
is a solution of the adjoint equation ^ to Eq. (7.35): 

L*[z]^0. (7.42) 

Knowledge of a solution (7.41) to the adjoint equation (7.42) allows to reduce the 
order of Eq. (7.35) by integrating Eq. (7.40): 

il:{x,y,y',...,y^^-'^) = C,. (7.43) 

Here Ci is an arbitrary constants and t/j defined according to Eqs. (7.31)-(7.32), 
i.e. 

ij = H{zL[y])l^^. (7.44) 

Proof. If (7.41) is a solution of the adjoint equation (7.42), we substitute it in 
Eq. (7.32) and arrive at Eq. (7.40). Hence (f){x) is an integrating factor for Eq. 
(7.35). Conversely, if (f){x) is an integrating factor for Eq. (7.35), then Eq. (7.40) 
is satisfied. Now Proposition 7.2 yields 

Si<P{x)L[y]) _ ^ 
5y 

Hence (7.41) is a solution of the adjoint equation (7.42). Finally, Eq. (7.44) 
follows from (7.32). 



^This statement is applicable to nonlinear ODEs as well, see [28]. 
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Example 7.2. Let us apply the above approach to the first-order equation 

y' + P{x)y = Q{x). (7.45) 
Here L[y] — y' + P{x)y. The adjoint equation (7.42) is written 

z' - P{x)z = 0. 
Solving it we obtain the integrating factor 

Therefore we rewrite Eq. (7.45) in the equivalent form 

[y' + P{x)y] e/^^^)"'^ = Q(a;)eJ'^(^)'^^ 
and compute the function ^ given by Eq. (7.31): 

* = U{zL[y\) = w-^\ziy' + P{x)y)\ ^wz^ «;eJ'^^^)'^^. 
dy' 



(7.46) 



Eq. (7.44) yields 



ij = ye 



f P{x)dx 



(7.47) 



Now we can take (7.46) instead of Eq. (7.38) and write it in the form (7.39) with 
ip given by (7.47). Then we obtain 



ye 



J P{x)dx 



j g(x)e/^(- 



= 0, 



whence 



y^fp{x)dx = Ci + y" g(x)e^^(^)'^^dx. 

Solving the latter equation for y we obtain the general solution of Eq. (7.45): 

/g(.,e/-<..-^; 



Ci + 



- / P{x)dx 



Example 7.3. Let us consider the second-order homogeneous equation 



sma; 



y + 



x^ 



COS X sm X 



x^ x^ 



y = 0. 



(7.48) 



(7.49) 



Its left-hand side does not satisfy the total derivative condition (7.19) because 



S_ 
5y 



y" + 



x^ 



y' + 



cos X sm X 



x^ 



x-' 



smx 



x^ 
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Therefore we will apply Proposition 7.4. The adjoint equation to Eq. (7.49) is 
written 

„ sin a; , sinx 
z" — z' + ^ = 0. 

X'^ X"^ 

We take its obvious solution z — x, substitute it in Eq. (7.31) and using (7.33) 
find 



^ = N 

Therefore Eq. (7.43) is written 

xy' + 



,, sm a; , / cos x sm x 

xy + y + Y- ) y 

X \ X X'^ 



X 



w — w + xw' . 



smx 

x 



Integrating this first-order linear equation we obtain the general solution of Eq. 
(7.49): 

C2 + / 1 e^ ^ ''^dx) xe- ^"^'^ (7.50) 



7.5 Application of the operator identity to linear PDEs 

Using the operator identity (7.1) one can easily extend the equations (7.32)- 
(7.33) for linear ODEs to linear partial differential equations and systems. Let 
us consider the second-order linear operator 



L[u\ — a^-'{x)uij + U'{x)ui -\- c{x)u 
considered in Section 1.5, Remark 1.1. The adjoint operator is 

L*[v\ = ^^^^^^ = DiDjia'^v) - D0v) + cv. 
Let us take the operator identity (7.1), 



(7.51) 



(7.52) 



X ^W— + DiN\ 
6u 

where X is the operator (7.14) with one dependent variable 

d d d 

X = W— + Wi— + W,^--, 
ou aui ouij 



and N* are the operators (7.5), 

— + Wj — 



Sui 



Suij 



W 



d_ 

dui 



ou 



d 



dui 



(7.53) 



52 



We use above the notation Wj = Di{W), Wij = DiDj{W). Now we proceed as 
in Section 7.3. Namely, we act on vL[u] by both sides of the identity (7.53), 

X{vL[u]) = W^-^^ + DiW{vL[u]), 
du 

take into account that X does not act on the variables x\ v, and that X{L[u]) = 
L[W], use Eq. (7.52) and obtain: 

vL[W]) - WL*[v] = DiN'{vL[u]). 

Letting here W — u we arrive at the following generalization of the equation 
(7.32): 

vL[u] - uL*[v] = Diiijj'), (7.54) 
where -0' are defined as in (7.33)-(7.34): 

ilj' = N\vL[u])\^^^ = a'^{x)[vui - uvi] + [b\x) - Di{a'^ {x))]uv. (7.55) 

7.6 Application of the operator identity to nonlinear equa- 
tions 

Let us apply the constructions of Section 7.5 to nonlinear equations (1.6), 

. . . = 0, Q; = l,...,m. (7.56) 

We write the operator (7.14) in the form 

Q Q O 

X = W"— + W"'— + W^- \ 

where P^" = Di{W"), W^^ = A^i(W^"), • • • • Then the operator (7.5) is written 

We act on v^Fp by both sides of the operator identity (7.1) 

X^W-^ + D.W, 

denote by F*[v] the adjoint operator defined by Eq. (1.8) and obtain 

v^Fp[W] - W^F^iv] = Di{¥), (7.57) 

53 



where 

and is the linear approximation to defined by (see also Section 1.3) 

r)P rlW rlW 

Fb\W] = X{Fb) = W—^ + W—^ + W"—^ + ■■■ . 

Remark 7.5. Eq. (7.57) shows that F*[v\ = i.e. the adjoint operator F* 

to nonlinear Eqs. (7.56) is the usual adjoint operator F^ to the linear operator 
(see also [29]). But the linear self-adjointness of -F/jfVT] is not identical 
with the nonlinear self-adjointness of Eqs. (7.56). For example, the KdV equation 

= is nonlinear ly self- adjoint (see Example 1.2 in Section 
1.6). But its linear approximation F[W] = Wf — W^xx ~ uW^ — Wu^ is not a 
self-adjoint linear operator. Moreover, all linear equations are nonlinearly self- 
adjoint. 



8 Conservation laws: Generalities and explicit 
formula 

8.1 Preliminaries 

Let us consider a system of m difi:erential equations 

= 0, Q; = l,...,m, (8.1) 

with m dependent variables m^, . . . , and n independent variables a;^, . . . , x". 
A conservation law for Eqs. (8.1) is written 

[A(C0](3.,) = O. (8.2) 

The subscript |(8 i) means that the left-hand side of (8.2) is restricted on the 
solutions of Eqs. (8.1). In practical calculations this restriction can be achieved 
by solving Eqs. (8.1) with respect to certain derivatives of u and eliminating these 
derivatives from the left-hand side of (8.2). For example, if (8.1) is an evolution 
equation 

Ut = ^(t,X,U,Ux,Uxx), 

the restriction |(8.i) can be understood as the elimination of Ut- The n-dimensional 
vector 

C = {C\...,C'') (8.3) 
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satisfying Eq. (8.2) is called a conserved vector for the system (8.1). If its 
components are functions C* = C"(a;, -u, ■U(i), . . .) oi x^u and derivatives ^(i), . . . 
of a finite order, the conserved vector (8.3) is called a local conserved vector. 

Since the conservation equation (8.2) is linear with respect to C\ any linear 
combination with constant coefficients of a finite number of conserved vectors is 
again a conserved vector. It is obvious that if the divergence of a vector (8.3) van- 
ishes identically, it is a conserved vector for any system of differential equations. 
This is a trivial conserved vectors for all differential equations. Another type of 
trivial conserved vectors for Eqs. (8.1) are provided by those vectors whose com- 
ponents C* vanish on the solutions of the system (8.1). One ignores both types 
of trivial conserved vectors. In other words, conserved vectors (8.3) are simplified 
by considering them up to addition of these trivial conserved vectors. 

The following less trivial operation with conserved vectors is particularly use- 
ful in practice. Let 

C'l(8.1) = + ^2(/^') + • • • + Dr,{H^) (8.4) 
the conserved vector (8.3) can be replaced with the equivalent conserved vector 

C= (C^C'^...,C") = (8.5) 

with the components 

C\ (5' = C' + Di(i/'), C"" = C'' + Di{H''). (8.6) 

The passage from (8.3) to the vector (8.5) is based on the commutativity of the 
total differentiations. Namely, we have 

and therefore the conservation equation (8.2) for the vector (8.3) is equivalent to 
the conservation equation 

'd,{&)\ = 

L '\ (8.1) 

for the vector (8.5). If n > 3, the simplification (8.6) of the conserved vector can 
be iterated: if C'^ contains the terms 

D^{H^) + ■ ■ ■ + D^{H^) 
one can subtract them from and add to C^, . . . , C" the corresponding terms 

D2(^'),...,i^2(^"). 
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Note that the conservation law (8.2) for Eqs. (8.1) can be written in the form 

D,{C')=ii"F^{x,u,U(^),...,U(^s)) (8.7) 

with undetermined coefficients fi" = /i"(x, u, ^(i), . . .) depending on a finite num- 
ber of variables x, u, ... . If C* depend on higher-order derivatives, Eq. (8.7) 
is replaced with 

A(C") = /x"F^ + //-A (Fa) + f^'^'DiDj (Fe,) + • • • . (8.8) 

It is manifest from Eq. (8.7) or Eq. (8.8) that the total differentiations of a 
conserved vector (8.3) provide again conserved vectors. Therefore, e.g. the vector 

Dr{C) = {D,{C'),...,D,{Cn) (8.9) 

obtained from a known vector (8.3) is not considered as a new conserved vector. 

If one of the independent variables is time, e.g. = t, then the conservation 
equation (8.2) is often written, using the divergence theorem, in the integral form 



/ dx^ • • • dx" = 0. (8.10) 



d 

dt JiR" 



But the differential form (8.2) of conservation laws carries, in general, more in- 
formation than the integral form (8.10). Using the integral form (8.10) one may 
even lose some nontrivial conservation laws. As an example, consider the two- 
dimensional Boussinesq equations 

Vt + fi>z = ll^xVz - i>zVx , (8.11) 

H V^x = H^xPz - PzPx 

9 

used in geophysical fluid dynamics for investigating uniformly stratified incom- 
pressible fiuid fiows in the ocean. Here A is the two-dimensional Laplacian, 

+ 



and is the stream function so that the x, z- components w of the velocity 
{u, V, w) of the fiuid are given by 

u = ipz, w = -ip^. (8.12) 

Eqs. (8.11) involve the physical constants: g is the gravitational acceleration, / 
is the Coriolis parameter, and is responsible for the density stratification of 



56 



the fluid. Each equation of the system (8.11) has the conservation form (8.2), 
namely 



Dt{AiP) + D,{~gp + ^,A^) + D,{-fv - ^,A^) = 0, 

Dt{v) + D,{vij,) + D,{f^P - v^P^) = 0, (8.13) 

A(P) + {— ,l; + p,l;A+ D,(-p^^) = 0. 

In the integral form (8.10) these conservation laws are written 

d 



J J A^j dxdz = 0, 



diJJ 



dxdz = 0, (8.14) 
pdxdz — 0. 



We can rewrite the differential conservation equations (8.13) in an equivalent 
form by using the operations (8.4)-(8.6) of the conserved vectors. Namely, let us 
apply these operations to the first equation (8.13), i.e. to the conserved vector 

C' = A,p, C^ = -gp + i^,A,p, C' = -fv-,p,A,p. (8.15) 

Noting that 

and using the operations (8.4)-(8.6) we transform the vector (8.15) to the form 

C' = 0, C^ = -gp + ^t, + ^,A^, d^ = -fv + iJu-i^xAij. (8.16) 

The integral conservation equation (8.10) for the vector for (8.16) is trivial, = 0. 
Thus, after the transformation of the conserved vector (8.15) to the equivalent 
form (8.16) we have lost the first integral conservation law in (8.14). But it does 
not mean that the conserved vector (8.16) has no physical significance. Indeed, 
if write the differential conservation equation with the vector (8.16), we again 
obtain the first equation of the system (8.11): 

D,{C^) + D,{C^) = Ai/Jt - gp, - fv, - V.A^, + ^,A^, . 

Let us assume that Eqs. (8.1) have a nontrivial local conserved vector sat- 
isfying Eq. (8.7). Then not all p^ vanish simultaneously due to non-triviality 
of the conserved vector. Furthermore, since p^Fg depends on x, u and a finite 
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number of derivatives ^(i), M(2), . . . (i.e. it is a differential function) and has a 
divergence form, the following equations hold (for a detailed discussion see [4], 
Section 8.4.1): 

(5 r " 

— /i''F^(a;,M,M(i), . . . =0, a = l,...,m. (8.17) 

Note that Eqs. (8.17) are identical with Eqs. (3.2) where the differential sub- 
stitution (3.9) is made with (p" = Hence, the system (8.1) is nonlinearly 
self-adjoint. 1 formulate this simple observation as a theorem since it is useful in 
applications (see Section 11). 

Theorem 8.1. Any system of differential equations (8.1) having a nontrivial 
local conserved vector satisfying Eq. (8.7) is nonlinearly self-adjoint. 



8.2 Explicit formula for conserved vectors 

Using Definition 3.1 of nonlinear self-adjointness and the theorem on conserva- 
tion laws proved in [3] by using the operator identity (7.1), we obtain the explicit 
formula for constructing conservation laws associated with symmetries of any 
nonlinearly self-adjoint system of equations. The method is applicable indepen- 
dently on the number of equations in the system and the number of dependent 
variables. The result is as follows. 

Theorem 8.2. Let the system of differential equations (8.1) be nonlinearly self- 
adjoint. Specifically, let the adjoint system (3.2) to (8.1) be satisfied for all 
solutions of Eqs. (8.1) upon a substitution (3.3), 



v° = ^'^{x, u), 



a 



1, . . . ,m. 



(8.18) 



Then any Lie point, contact or Lie-Backlund symmetry 



d d 



^.19) 



as well as a nonlocal symmetry of Eqs. (8.1) leads to a conservation law (8.2) 
constructed by the following formula: 



dC 



dC 



du 



+ D.Dk 



dC 



(8.20) 



dC 



dC 



+ . . . 



+ DjDk 



dC 
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where 

and jC is the formal Lagmngian for the system (8.1), 



(8.21) 



(8.22) 



In (8.20) the formal Lagrangian C should be written in the symmetric form with 
respect to all mixed derivatives m^-, m^-^,, . . . and the "non- physical variables" v" 
should be eliminated via Eqs. (8.18). 

One can omit in (8.20) the term when it is convenient. This term provides 
a trivial conserved vector mentioned in Section 8.1 because C vanishes on the 
solutions of Eqs. (8.1). Thus, the conserved vector (8.20) can be taken in the 
following form: 



dC 

du? 



dC 



dC 



du 



ijk 



(8.23) 



+ D, (W") 



dC 



dC 



+ . . . 



+ DjDk {Wn 



dC 



du 



ijk 



Remark 8.1. One can use Eqs. (8.23) for constructing conserved vectors even 
if the system (8.1) is not self-adjoint, in particular, if one cannot find explicit 
formulae (8.18) or (3.9) for point or differential substitutions, respectively. The 
resulting conserved vectors will be nonlocal in the sense that they involve the 
variables v connected with the physical variables u via differential equations, 
namely, adjoint equations to (8.1). 

Remark 8.2. Theorem 8.2, unlike Nother's theorem 7.1, does not require addi- 
tional restrictions such as the invariance condition (7.9) or the divergence condi- 
tion mentioned in Remark 7.2. 



9 A nonlinearly self-adjoint irrigation system 

Let us apply Theorem 8.2 to Eq. (6.1) satisfying the condition (6.3): 

C(V')^t = [i^(^)^.], + im) (^. - 1)], - ^(V^), (9.1) 
S\ip) = aC{ip), a = const. (9.2) 

The formal Lagrangian (8.22) for Eq. (9.1) has the form 

C = [-C{il^)il^t + i^(^)(^.. + ^,,) + i^'(^)(^^ + il^l - ^.) - 5(^)] V. (9.3) 
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We will use the substitution (6.5) of the particular form 



(9.4) 



Denoting t = x^,x — x^,z — we write the conservation equation (8.2) in the 
form 

A(C^) + D,{C^) + D,{C') = 0. (9.5) 

This equation should be satisfied on the solutions of Eq. (9.1). 

The formal Lagrangian (9.3) does not contain derivatives of order higher than 
two. Therefore in our case Eqs. (8.23) take the simple form 



C' = W 



dC 



dC 



(9.6) 



and yield: 



dC 

dipt 



W 
W 









dipx 




yd'ipxxj 


' dC 




( dc Y 






dipz 





+ D,{W) 



dC 

dC 



Substituting here the expression (9.3) for C we obtain 

= -WC{ij)v, 

= W[2K\i,)vij, - D^{K{i,)v)] + D^{W)K{i,)v, 

= W[K\i,)v{2ilj, - 1) - D.,{K{^j)v)] + D.XW)K{ij)v, 

where v should be eliminated by means of the substitution (9.4). So, we have: 

= -WC{'iP)e''\ 

C = [WK\i^)i^, + D,{W)Kme''\ (9.7) 
= [WK\iP){iP, - 1) + D,{W)K{iP)]e^\ 

Since Eq. (9.1) does not explicitly involve the independent variables t, x, it is 
invariant under the translations of these variables. Let us construct the conserved 
vector (9.7) corresponding to the time translation group with the generator 



X- ^ 



(9.8) 



For this operator Eq. (8.21) yields 



(9.9) 
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Substituting (9.9) in Eqs. (9.7) we obtain 

C2 = -[ir'(V)Vt^. + K{^),Pt.]e''\ (9.10) 
= -[K'{ilj)ilJt{ilj, - 1) + i^(^)V^t,]e«*. 

Now we replace in the term C{'^)'^t by the right-hand side of Eq. (9.1) to 
obtain: 

= -5(^)e"* + (i^(^)^.e'^*) + {K{^P){^P, - l)e"*) . 

When we substitute this expression in the conservation equation (9.5), we can 
write 

Therefore we can transfer the terms Dx(. . .) and Dz{- ■ ■) from to and C^, 
respectively (see (8.6)). Thus, we rewrite the vector (9.10), changing its sign, as 
follows: 

= -S(V')e"*, 

= [ir'(V')V'tV'x + i^(V')V'tx]e"* - A (i^(V')V'xe"*) , 

= [K'{i^)i^t{i^, - 1) + X(V')V't.]e"* - A (i^(V')(V'. - l)e"*) . 

Working out the differentiation Dt in the last terms of and we finally arrive 
at the following vector: 

C = aK(V')Vxe"*, (9.11) 
C^ = aX(V')(V'.-l)e"*. 

The reckoning shows that the vector (9.11) satisfies the conservation equation 
(9.5) due to the condition (9.2). Note that is the density of the conserved 
vector (9.11). 

The use of the general substitution (6.5) instead of its particular case (9.4) 
leads to the conserved vector with the density 

= 5(^/;)(6a; + 0e"*. 

This approach opens a new possibility to find a variety of conservation laws 
for the irrigation model (6.1) by considering other self-adjoint cases of the model 
and using the extensions of symmetry Lie algebras (see [14], vol. 2, Section 9.8). 
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10 utilization of differential substitutions 



10.1 Equation Uxy = sinu 

We return to Section 3.2 and calculate the conservation laws for Eq. (3.11), 

Uxy — sinii, (10-1) 
using the differential substitution (3.12), 

V = Ax\xux - yu^ + AiUx + A-^Uy, (10-2) 

and the admitted three-dimensional Lie algebra with the basis 

d (9 B B 

Xi^^, X2 = —, X^^x—-y—- (10.3) 

ox oy ox oy 

The conservation equation for Eq. (10.1) will be written in the form 

DxiC') + Dy{C^)^0. 
We write the formal Lagrangian for Eq. (10.1) in the symmetric form 

jO-= Q Mxj/ + ^ Uyx - sin it^ V. (10.4) 



Eqs. (8.23) yield: 
1 



= - Dy{W)v - - Wvy, C^ = - Dx{W)v - - Wvx. (10.5) 



where we have to eliminate the variable v via the differential substitution (10.2). 

Substituting in (10.5) W = —u^ corresponding to the operator Xi from (10.3), 
replacing v with (10.2) and u,j.y with sin-u, then transferring the terms of the form 
Dy{. . .) from to (see the simplification (8.6)) we obtain: 

= ^licosti, C'^^l-Aiul- 

We let Ai = 1 and conclude that the application of Theorem 8.2 to the symmetry 
Xi yields the conserved vector 

C^^cosu, C'^^^ul- (10.6) 

The similar calculations with the operator X2 from (10.3) lead to the conserved 
vector 

C^ = ^ul, C^^cosu. (10.7) 
62 



The third symmetry, X3 from (10.3), does not lead to a new conserved vector. 
Indeed, in this case W = yUy — xu^- Substituting it in the first formula (10.5) we 
obtain after simple calculations 

= ^Asul - A2 COSU + Dy 

Hence, upon transferring the term Dy{. . .) from to the resulting will 
be a linear combination with constant coefficients of the components of the 
conserved vectors (10.6) and (10.7). The same will be true for . Therefore the 
conserved vector provided by the symmetry X3 will be a linear combination with 
constant coefficients of the conserved vectors (10.6) and (10.7). 

One can also use the Noether theorem because Eq. (10.1) has the classical 
Lagrangian, namely 

L — —-Uy:Uy^ Q,o%u. (10.8) 

Then the symmetries X\ and X2 provide again the conserved vectors (10.6) and 
(10.7), respectively. But now we obtain one more conserved vector using X3, 
namely 

= xcosii- C"^ ^^ul-yco%u. (10.9) 
10.2 Short pulse equation 

The differential equation (up to notation and appropriate scaling the physical 
variables) 

DtD,{u)^u^^-Dl{u'-) (10.10) 

was suggested in [30] (see there Eq. (11), also [31]) as a mathematical model 
for the propagation of ultra-short light pulses in media with nonlinearities, e.g. 
in silica fibers. The mathematical model is derived in [30] by considering the 
propagation of linearly polarized light in a one-dimensional medium and assuming 
that the light propagates in the infrared range. The final step in construction of 
the model is based on the method of multiple scales. 

Eq. (10.10) is connected with Eq. (10.1) by a non-point transformation 
which is constructed in [32] as a chain of differential substitutions (given also in 
[31] by Eqs. (2)). Using this connection, an exact solitary wave solution (a pulse 
solution) to Eq. (10.10) is constructed in [31]. One can also find in [32] a Lax 
pair and a recursion operator for Eq. (10.10). 

Note that Eq. (10.10) does not have a conservation form. I will find a conser- 
vation law of Eq. (10.10) thus showing that it can be rewritten in a conservation 
form. A significance of this possibility is commonly known and is not discussed 
here. 



{A2y + A^x) ( -UxUy + cosM 
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We write the short pulse equation (10.10) in the expanded form 

Uxt^u+^ u^Uxx + uul 
so that the formal Lagrangian is written 



Uxt ~ ~ 2 ^ 



:io.ii) 



(10.12) 



Substituting (10.12) in (3.2) we obtain the following adjoint equation to Eq. 
(10.11): 

1 



Vxt=V + - u Vx 



(10.13) 



We first demonstrate the following statement. 
Proposition 10.1. Eq. (10.10) is not nonlinearly self-adjoint with a substitution 

v = ip(t,x,u) (10.14) 

but it is nonlinearly self- adjoint with the differential substitution 

1 . 



V = Ut - -U Ux- 



;i0.15) 



Proof. We write the nonlinear self-adjointness condition (3.5), 



Vxt-V- - U Vxx 



(10.14) 



X[uxt - - 2 ^^^^^ ~ 



substitute here the expression (10.14) for v and its derivatives 

'^xx — ^ufJ'xx ~l~ ^uufJ'x ~^ '^Vxu'^x ~l~ ^XXl 

Vxt = <PuUxt + <PuuUxUt + (PxuUt + <PtuUx + <Pxt , 



;io.i6) 



and first obtain A = by comparing the terms with the second-order deriva- 
tives of u. This reduces the nonlinear self-adjointness condition to the following 
equation: 



^uuUxUt + ^xuUt + ^tuUx + ^xt-^- T; u'^i^'uuUl + ^if^^Ux + ^xx) 



(10.17) 



= -iPu[u + Uul]. 

The terms with Ut in Eq. (10.17) yield 99 
with ul and obtain = 0. Hence 



UU TXU 



0. Then we take the term 



ip = a{t, x). 
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Now Eq. (10.17) gives = 0, a^t — a — 0, whence a = 0. Thus 

(^ = 0, 

i.e. the substitution (10.14) is trivial. This proves the first part of Proposition 
10.1. Its second part is proved by similar calculations with the substitution 

V = (p{t,X,U,Ux,Ut). 

I will not reproduce these rather lengthy calculations, but instead we will verify 
that the substitution (10.15) maps any solution of Eq. (10.1) into a solution of 
the adjoint equation (10.13). First we calculate 

Vx =Uxt-^ U^Uxx - Uul 

and see that on the solutions of Eq. (10.1) we have Vx — u. Now we calculate 
other derivatives and verify that on the solutions of Eq. (10.1) the following 
equations hold: 

Vx ^U, Vt^ Utt - ^ U^Uxt - UUxUt, Vxt = Ut, Vxx = Ux- (10.18) 

It is easily seen that Eq. (10.13) is satisfied. Namely, using (10.15) and (10.18) 
we have: 

Vxt — V U Vxx = Uf— \ Ut U Ux \ U Ux = 

The maximal Lie algebra of point symmetries of Eq. (10.10) is the three- 
dimensional algebra spanned by the operators 

d d d d d 

at ox ou ox at 

Let us construct the conservation laws 

A(C') + ^x(O = (10.20) 

for the basis operators (10.19). 

Since the formal Lagrangian (10.12) does not contain derivatives of order 
higher than two, Eqs. (8.23) are written (see (9.6)) 



dui ^ \duij 
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+ DAW) 



duij 



In our case we have: 
C' = -WDj 

dC 



dU-r 



dutx 
-A 



dC 



du 



tx 



dC 



du 



xt 



dC 



dC 



(10.21) 



Substituting in (10.21) the expression (10.12) for £ written in the symmetric form 

11 1 o 



^Utx + Uxt 



U — -U Ua 

2 



we obtain 



(10.22) 



(10.23) 



= -W 



1 



1 



UVUx H — Vt U Vx 

2 2 



Since v should be eliminated via the differential substitution (10.15), we further 
simplify this vector by replacing Vx with u according to the first equation (10.18) 
and obtain: 



C^^-^^Wu+^^vDxiW), 



(10.24) 



= -W 



uvux + ^ ^* - ^ ^""J + ^ vDt{W) - ^ u^DxiW), 



where v and Vt should be replaced with their values given in Eqs. (10.15), (10.18). 

Let us construct the conserved vectors using the symmetries (10.19). Their 
commutators are 

[Xi.Xs] = —Xi, [X2,X3] = X2. 

Hence, according to [10], Section 22.4, the operator X3 plays a distinguished role. 
Namely, the conserved vectors associated with Xi and X2 can be obtained from 
the conserved vector provided by X3 using the adjoint actions of the operators 
Xi and X2, respectively. Therefore we start with X3. Substituting in (10.24) the 
expression 

W = U + tUt — XUx 

corresponding to the symmetry X3, eliminating the terms of the form D,j.{A) from 
and adding them to in the form Dt{A) according to the simplification (8.6), 
we obtain after routine calculations the following conserved vector: 



u 



2 2-^4-^42 
U UxUt — Uf — -U ~ '^x- 



(10.25) 
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The conservation equation (10.20) for the vector (10.28) holds in the form 

A(C') + D^{C^) = 2(^ut-^ u\^^ + i u\^^ + uul - u^t) ■ (10.26) 

Let us turn now to the operators Xi and X2 from (10.19). To simplify the 
calculations it is useful to modify Eqs. (10.24) as follows. Noting that 

vD,{W) = D,{vW) - Wv, 

we rewrite the vector (10.23) in the form 

= -W 



2 



+ vDt(W)-]^u'vD,{W). 



Then (10.24) is replaced with 

C-i = -uW, (10.27) 

= -W 



1 3 

UVUj; — - U 



+ vDt{W)-^u''vD,{W). 



Substituting in the first formula (10.27) to expression W — —Ut corresponding the 
operator Xi we obtain — uUf. This is the time derivative of from (10.25). 
Hence the symmetry Xi leads to a trivial conserved vector obtained from the 
vector (10.25) by the differentiation Dt, in accordance with [10]. Likewise, it is 
manifest from (10.27) that the operator X2 leads to a trivial conserved vector 
obtained from the conserved vector (10.25) by the differentiation D^.. Thus we 
have demonstrated the following statement. 

Proposition 10.2. The Lie point symmetries (10.19) of Eq. (10.11) yield one 
non-trivial conserved vector (10.25). Accordingly, the short pulse equation (10.11) 
can be written in the following conservation form: 

Dt (u'^^ + (u^u^ut - ul - ^u" - ^M^) = 0. (10.28) 

11 Gas dynamics 

11.1 Classical symmetries and conservation laws 

Let us consider the polytropic gasdynamic equations 

Vt + (v • V)v + - Vp = 0, 
P 

Pt + v-Vp + pV -v^O, (11.1) 
Pt + v ■ Vp + 7pV -v — O, 
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where 7 is a constant known as the polytropic (or adiabatic) exponent. The 
independent variables are the time and the space coordinates: 

t, a; = n<3. (11.2) 

The dependent variables are the velocity, the density and the pressure: 

v = {v\...,v''), p, p. (11.3) 

Eqs. (11.1) with arbitrary 7 have the Lie algebra of point symmetries spanned 

by 

X — ^ X — ^ Y —t^ Y- — t ^ ^ 

° 9i ' ' dx'' ' ° dt dx^ ' ' dx"^ dv'^ ' 

d d ^ d , d d 
^o^pTT+P^^ ^1 = -"t^ ^ + 2p^, z,j = l,...,n, 
op op ot ov^ op 

and the following classical conservation laws: 
d f 

— / pdcu — - Conservation of mass 



dt 

n{t) 

d 



j (^^''^^^^^^ ^^—^duj — — Jpv-udS - Energy 



dt I l^^l" 
n(t) 5(t) 

d 



, , pvduj = — p fdS - Momentum 

dt J J 

n{t) s{t) 

— J p{x X v)du! — — J p{x X v)dS - Angular momentum 

^ j p{tv — x)duj — — J tpvdS - Center-of-mass theorem. 

s{t) 

The conservation laws are written in the integral form by using the standard 
symbols: 

Vt{t) - arbitrary n-dimensional volume, moving with fluid, 

S{t) - boundary of the volume n(t), 

v - unit (outer) normal vector to the surface S{t). 

If we write the above conservation laws in the general form 

^[Tduj=-Ux- ^)dS, (11.5) 



dt „ 

ci{t) sit) 
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then the differential form of these conservation laws will be 

Dt{T) + V-{x + Tv) = 0. 



(11.6) 



11.2 Adjoint equations and self-adjointness when n = 1 

Theorem 8.1 from Section 8.1 shows that the system of gasdynamic equations 
(11.1) is nonlinearly self-adjoint. Let us illustrate this statement in the one- 
dimensional case: 

Vt + vv^ + = 0, 
P 

Pt + Vpa:+ pv,:^0, (11.7) 
Pt + VPx + Wt^x = 0. 

We write the formal Lagrangian in the form 

C = u(vt + vvx + -^Pxj+ ^(Pt + vpx + pVx) + P{Pt + vpx + jpvx) (11.8) 
and obtain the following adjoint system for the new dependent variables U,R,P: 
— = -Ut- vUx - pRx + (1 - 7)Ppx - IpPx = 0, 



5C ^ ^ 1 

— = -Rt - vR^ ; 

dp p- 

5^ ^ 1 1 
= -P^--Ux + -, 
dp p p- 

Let us take, e.g. the conservation of energy from Section 8.1. Then we have 



= -Rt -vR^-- Up^ = 0, (11.9) 
'^ = -Pt-tu, + \upx + i^- l)Pvx - vPx = 0. 



T^-pv H , X = pv, 

2 7 — 1 

and using the differential form (11.6) of the energy conservation we obtain the 
following equation (8.7): 

= PV (vt + VVx + -Px)+^iPt + Vpx + PVx)-^ ^ {Pt + VPx + IPVx)- 

\ p J I 7 — 1 

Hence, the adjoint equations (11.9) are satisfied for all solutions of the gasdynamic 
equations (11-1) upon the substitution 

U^pv, R^%. P-^- (11-11) 
z 7 — i 
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This conclusion can be easily verified by the direct substitution of (11.11) in the 
adjoint system (11.9). Namely, we have: 



5C 

5v 

SC 

Tp 

SC 
Sp 



(11.11) 



(11.11) 



= -p(vt + vv^c ^~P^)~ '"(Pt + '"P=' + P^^)^ 

-Px), 
o / 



(11.11) ^ p 

= -vivt + Wfa; + 



0. 



;ii.i2) 



11.3 Adjoint system to equations (11.1) with n > 2 

For gasdynamic equations (11.1) with two and three space variables the formal 
Lagrangian (11.8) is replaced by 



C^U ■ (vt + {v- V)v + ^ Vp) + R{pt + -y • Vp + pV • v) 

+ P{pt + v - Vp + ^^pV -v), (11.13) 



where the vector U = (f/^, . . . , f/") and the scalars R,P are new dependent 
variables. Using this formal Lagrangian, we obtain the following adjoint system 
instead of (11.9): 



— =-Ut-{v- V)U + {U ■ V)v - (V • v)U 



Sv 



- pVR + (1 - 7)PVp - 7pVP = 0, 



— =-Rt-v-VR- — U-Vp = Q, (11.14) 
op p^ 

Ft r 1 1 

— = - Pi - - (V •[/) + -[/• Vp+ (7 - 1)P(V • v) - V • VP = 0. 

bp p p^ 



The nonlinear self-adjointness of the system (11.1) can be demonstrated as in the 
one-dimensional case discussed in Section 11.2. 
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11.4 Application to nonlocal symmetries of the Chaplygin 
gas 

The Chaplygin gas is described by the one-dimensional gasdynamic equations 
(11.7) with 7 = -1 : 

vt + vvx + -p^ = 0, 
P 

Pt + Vpj:+ pV^^O, (11.15) 
Pt + VPx - PVx = 0. 

Eqs. (11.15) have the same maximal Lie algebra of Lie point symmetries as Eqs. 
(11.7) with arbitrary 7. This algebra is spanned by the symmetries (11.4) in the 
one-dimensional case, namely 

d d d d d d 

at ox ox ov ot ox 

d d d d d , ^ 

But the Chaplygin gas has more symmetries than an arbitrary one-dimensional 
polytropic gas upon rewriting it in Lagrange's variables obtained by replacing x 
and p with y and q, respectively, obtained by the following nonlocal transforma- 
tion: 

T= [ pdx, - • (11-17) 

J P 

Then the system (11.15) becomes 

qt-Vr^ 0, 

Vt+Pr^O, (11.18) 

Pt - -Vr^O 
Q 

and admits the 8-dimensional Lie algebra with the basis 

9 9 9 (9 d 

ot OT OV ot OT 

d d d d d d 

n = r— +p— , Ye^v—+p— + q—, 11.19 
OT Op Oq Ov Op Oq 

d q d d d yq d 

^7 = ^ + - ^ > ^8 ^ ~ ^ ■ 

op p oq ov op p oq 
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It is shown in [21] that the operators Y-^,Ys, from (11.19) lead to the following 
nonlocal symmetries for Eqs. (11.15): 



^ d d p d 

^ dx dp p dp^ 

X-l- s^— t—-T— ^ — 
^ 2 J dx dv dp p dp^ 



;il.20) 



where T,s,a are the following nonlocal variables: 

T= [ pdx, s = — I — dx, cr = — [ — ■ (11.21) 
J J P J P 

They can be equivalently defined by the compatible over-determined systems 

Tx = P, Tt+ VTa: = 0, 

Sx = --, St + VSx = 0, (11.22) 
P 

(Tx = , Crt + V(Tx = 0, 



p 



or 





T 


Tt = ■ 


-vp, 

VT 




P 


St = 








P 




1 




V 




P 










P 



(11.23) 



Let us verify that the operator is admitted by Eqs. (11.15). Its first prolon- 
gation is obtained by applying the usual prolongation procedure and eliminating 
the partial derivatives a^. and at via Eqs. (11.23). It has the form 

^ ^ ^ + ^ _^ _ _^ + _^ t^Px d ^ px d 

^ dx dp p dp p dvt p dvx p dpt p dp^ 

\p p^ P J dpt \ p p^ J dpx 

The calculation shows that the invariance condition is satisfied in the following 
form: 

^7 (vt + VV^+^ Px^ = 0, 

Xript + vp^ + pv^) = - (pt vp^ pv^) - ^{pt + vp^ - pv^), 

p p^ 

X'jipt + vpx-pvx) = 0. 
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One can verify likewise that the invariance test for the operator Xg is satisfied in 
the following form: 



Xg (^t + vvx + ^Px^ = 0, 

T OT 
Xs{pt + Vpa; + pVx) ^ -{Pt + VPx + pVx) 5- (Pt + VPx " PVx), 

XsiPt + vpx - pVx) = 0. 

The operators Yi, . . . ,Yq from (11.19) do not add to the operators (11.16) new 
symmetries of the system (11.15). 

Thus, the Chaplygin gas described by Eqs. (11.15) admits the eight-dimensional 
vector space spanned by the operators (11.16) and (11.20). However this vector 
space is not a Lie algebra. Namely, the commutators of the dilation generators 
X4, X5, Xq from (11.16) with the operators (11.20) arc not linear combinations 
of the operators (11.16), (11.20) with constants coefficients. The reason is that 
the operators X4,X5,Xq are not admitted by the differential equations (11.22) 
for the nonlocal variables T,s,a. Therefore I will extend the action of the dila- 
tion generators to r, s, a so that the extended operators will be admitted by Eqs. 
(11.22). 

Let us take the operator X^. We write it in the extended form 

,^,9 d d ^ d d 
X'4 = t— + X— + a— + f3— + ^i— , 
at ox OT OS oa 

where a, ^, p are unknown functions of x, v, r, s, a. Then we make the pro- 
longation of X4 to the first-order partial derivatives of the nonlocal variables with 
respect to t and x by treating T,s,a as new dependent variables and obtain 

,,,9 d d ^ d d 
X'4 = t— + x— + a— + /3— + p— 
Ot ox OT OS Oa 

+ [Dt{a) - ^*]^ + [Dx{(^) - ^a.]^ 
+ [D,{f3)-st]— + [DM-Sx]g^^ 
+ [Dt{p) -<Tt]^ + [DM - Ox 



o ' L / X I o 

oat oax 
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Now we require the invariance of Eqs. (11.22): 



X'^ \^s^ + - J = 0, X'^ist + vs^) = 0, 
Xi (^a^ + ^) = °' + v^x) = 0. 



;il.25) 



As usual, Eqs. (11.25) should be satisfied on the solutions of Eqs. (11.22). Let 
us solve the equations X'4(tx: — p) — 0, X'^{Tt + vt^) — 0. They are written 

[DM) - r.\ii.22) - 0, [Dt{a) -n + v {D,{a) - t,)\^^^^^-^ = 0. (11.26) 

Since — Dy.{oi), the first equation in (11.26) is satisfied if we take 

a — T 

With this a the second equation in (11.26) is also satisfied because Tj + vt^ — 0. 
Now the first equation in the second line of Eqs. (11.25) becomes 



DM - 



P 



DM - 2s, = 



and yields 



(11.22) 



= 2s. 



The second equation in the second fine of Eqs. (11.25) is also satisfied with this 
p. Applying the same approach to the third line of Eqs. (11.25) we obtain 



p = a. 

After similar calculations with X^ and Xq we obtain the following extensions 
of the dilation generators: 



d d d ^ d d 

dt dx dr ds da ' 

d d d d 

dp ^ dp ^ dr ^ da ' 



(11.27) 



d d ^ d ^ d ^ d 
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The operators (11.20), (11.27) together with the operators Xi, X-^ from 
(11.16) span the eight-dimensional Lie algebra Lg admitted by Eqs. (11.15) and 
Eqs. (11.22). The algebra Lg has the following commutator table: 





Xi 


X2 


^3 


X', 


X', 


X', 


Xj 


Xs 


Xi 








X2 


Xi 





Xi 





X3 


Xo 











Xo 
















-X2 














-Xs 








X', 


-Xi 


-X2 

















Xs 


X', 




















-Xr 





X'e 


-Xi 





Xs 














2X8 


Xj 














Xj 











Xs 


-Xs 








-Xs 





-2X8 









Let us apply Theorem 8.2 to the nonlocal symmetries (11.20) of the Chaplygin 
gas. The formal Lagrangian (11.8) for Eqs. (11.15) has the form 

jC^u(vt + vVx + ^ + R{pt + vpx + pvr,) + P(pt + vpr, - pvx). (11.28) 

Accordingly, the adjoint system (11.9) for the Chaplygin gas is written 
5C 

— = -Ut- vU^ - pR^ + 2Pp^ + pP^ - 0, 
dv 

^ = -Pi - vR, -\up, = 0, (11.29) 
dp p^ 

SjC 1 1 

= -P^--U,+ Up,- 2Pv, - vP, = 0. 
dp p p^ 

Let us proceed as in Section 11.2. Namely, let us first construct solutions to the 
adjoint system (11.29) by using the known conservation laws given in Section 8.1. 

Since the one-dimensional does not have the conservation of angular momentum, 
we use the conservation of mass, energy, momentum and center-of-mass and ob- 
tain the respective differential conservation equations (see the derivation of Eq. 
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(11.10)): 

Dt{p) + D^{pv) = pt + vpx + pv^, (11.30) 

Dt{pv'^ - p) + D^{pv + pv^) = 2pv (vt + VV^ + ^ Pa;) 

+ v'^{pt + vp^ + pv^) - {pt + vpx - pvx), (11.31) 

Dt{pv) + D^{p + pv'^) = p(vt + vv^ + ^p^^ 

+ v{pt + vpx + pv,:), (11.32) 
Dtitpv — xp) + Dx{tp + tpv^ — xpv) 

= tp(vt + vv^ + ^Px^ + {tv - x){pt + vpx + pv^). (11.33) 

Eqs. (11.30)- (11.33) give the following solutions to the adjoint equations (11.29): 

C/ = 0, i?=l, P = 0, (11.34) 

U^2pv, = P = -1, (11.35) 

C/ = p, = P = 0, (11.36) 

C/ = ip, R^tv-x, P = 0. (11.37) 

The formal Lagrangian (11.28) contains the derivatives only of the first order. 
Therefore Eqs. (8.23) for calculating the conserved vectors take the simple form 

We denote 

t — x^, X — x'^, V — u^, p — u^, p — u^. 
In this notation conservation equation (8.2) will be written in the form 

[A(C^) + i^.(C^')](n.i5) = 0- (11-39) 
Writing (11.38) in the form 

dvt opt opt 



9^;^; dp:c dp:c 
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and substituting the expression (11.28) for £ we obtain the following final ex- 
pressions for computing the components of conserved vectors: 

^UW^ + RW^ + PW^, (11.40) 

= {vU + pR- pP)W^ + vRW'^ + Q [/ + vP^ W^, (11-41) 

where 

Vr°^ = r/"-f<, a = 1,2,3. (11.42) 

We will apply Eqs. (11.40)-(11.41) to the nonlocal symmetries (11.20). First 
we write the expressions (11.42) for the operator X-j from (11.20): 

W^^-av,, W^^^-ap,, W^^-{l + ap,). (11.43) 

P 

Then we substitute (11.43) in (11.40)-(11.41) and obtain four conserved vectors 
by replacing U,R,P with each of four different solutions (11.34)-(11.37) of the 
adjoint system (11.29). Some of these conserved vectors may be trivial. We select 
only the nontrivial ones. 

Let us calculate the conserved vector obtained by eliminating U, R, P by using 
the solution (11.34), t/ = 0, = 1, P = 0. In this case (11.40)-(11.41) and 
(11.43) yield 

C'^W^^-ap,, (11.44) 
P 

= pW^ + vW"^ = -apva: + -V- avp:c- 

p 

We write 

-crp^ = -D^{ap) + pcTj;, 
replace with —1/p according to Eqs. (11.22) and obtain 

= -a,{ap). 

Therefore application of the operations (8.4)-(8.6) yields — Q and 

= -(jpv^ + -V- avpx - Dt{ap) 
P 

p 

= -apv^ + - V - avpx - apt - atp 
P 

= -(^{Pt + vpa: + pv:,). 

We have replaced at with v/p according to Eqs. (11.23). The above expression 
for vanishes on Eqs. (11.15). Hence, the conserved vector (11.44) is trivial. 
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Utilization of the solutions (11.35) and (11.36) also leads to trivial conserved 
vectors only. Finally, using the solution (11.37), 

U = tp, R^tv-x, P = 0, 

we obtain, upon simplifying by using the operations (8.4)-(8.6), the following 
nontrivial conserved vector: 

= ap, = apv + t. (11.45) 

The conservation equation (11.39) is satisfied in the following form: 

A(C^) + D,{C^) = a{pt + vp, + pv,). (11.46) 

Note that we can write in (11.45) without t since it adds only the trivial 
conserved vector with the components = 0, = t. Thus, removing t in 
(11.45) and using the definition of a given in (11.21) we formulate the result. 

Proposition 11.1. The nonlocal symmetry of the Chaplygin gas gives the 
following nonlocal conserved vector: 

C^ = -p[^^ C^^-pvf — . (11.47) 
J P J P 

Mow we use the operator Xg from (11.20). In this case 

^ + P., (11-48) 
- + sj p.. 

Substituting in (11.40)-(11.41) the expressions (11.48) and the solution (11.34) 
of the adjoint system, i.e. letting U = 0, R = 1, P = 0, we obtain 

C' = pW' + vW'=tp+^- (^ + s) (pf, + fPx). 





= t - 




- PL 




p 




— —T 



Noting that 



-[\ + s\p^^-P^-D^{'-p + ps 
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we reduce the above vector to the trivial conserved vector = 0, — 0. 
Taking the solution (11.35) of the adjoint system, i.e. letting 



U = 2fyv, R^v^, P = -1, 



we obtain 



P \2 



= t(3pv'^ +p) + - VT 
P 



+ s {^pv^v^ + v^px + pv^ + vp^). 



+ 5 ) {pv'^ - P) 



Then, upon rewriting in the form 

= 2tpv + 2r-D^ 

and applying the operations (8.4)-(8.6) we arrive at the following conserved vec- 
tor: 

C^ = tpv + T, C'^ = t{pv^ +p). (11.49) 
The conservation equation (11.39) is satisfied for (11.49) in the following form: 

A(C^) + A(C^) = tp (^vt + vv^ + ^ p^^ + tv{pt + vp^ + pv^). (11.50) 

Taking the solution (11.36) of the adjoint system, i.e. letting 

U = p, R = v, P = 0, 

we obtain 

= pW^ + vW^, = 2pvW^ + v^W^ + W\ 

Substituting the expressions (11.48) for VT^ VT^, VT"^ and simplifying as in the 
previous case we obtain the conserved vector 

= tp, = tpv - T. (11.51) 

The conservation equation (11.39) is satisfied for (11.49) in the following form: 

Dt{C') + D,{C^) = t{pt + vp, + pv,). (11.52) 
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Finally, we take the solution (11.37), U — tp, R — tv — x, P — 0, and obtain 
= tpW^ + {tv -x)W'^, C^{2tpv - xp)W^ + {tv^ - xv)W^ + tW^. 
Simplifying as above, we arrive at the conserved vector 



C'^\^--sjp, C-'=\^--sjpv-tr. (11.53) 

The conservation equation (11.39) is satisfied for (11.49) in the following form: 

A(C^) + D,{C^) = - iPt + + pv.). (11.54) 

Substituting in the conserved vectors (11.49), (11.51) and (11.53) the defini- 
tion (11.21) of the nonlocal variables we formulate the result. 

Proposition 11.2. The nonlocal symmetry Xs of the Chaplygin gas gives the 
following nonlocal conserved vectors: 



= tpv + J pdx, = t{pv'^ + p); 

= tpv — / pdx; 





= tp, 




















p 


















p 



pdx j dx 
pdx ) dx 



(11.55) 
(11.56) 

(11.57) 



pv — t pdx 



Theorem 11.1. Application of Theorem 8.2 to two nonlocal symmetries (11.20) 
gives four nonlocal conservation laws (11.47), (11.55)-(11.57) for the Chaplygin 
gas (11.15). 



11.5 The operator identity for nonlocal symmetries 

Example 11.1. Let us verify that the operator identity (7.1) is satisfied for the 
nonlocal symmetry of the Chaplygin gas. Specifically, let us check that the 
coefficients of 

ddddddddd . 
dv ^ dp^ 5p ' dvt ' dvx ' dpt ' dp^ ' dpt ' dp^ 
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in both sides of (7.1) are equal. Using the first prolongation (11.24) of Xy and 
the definition of the nonlocal variable a given in Eqs. (11.23) we see that the 
left-hand side of the identity (7.1) is written 



X + D (C) = a— ^ \ P ^ vvg, d ^ d vpx d 
^ ' dx dp p dp p dvt p dvx p dpt 



Px d 

+ — h 

P dpx 



Pt PPt vpx 



p p^ p J dpt 



d 



+ 2 



P 



p2 



d 



dp: 



1 

X P 



;ii.59) 



Then we use the expressions (11.43) of VF" for the operator X^, substitute them 
in the definition (7.5) of and obtain in our approximation: 



Ml 9 

dvt 



p \ d , (9 

P J dpt 



a — av^ 



_d_ 
dv^ 



+ { - - (^P3 
P 



_d_ 
dpx 



dpt ' 
dpt 



Now the right-hand side of (7.1) is written: 



dv 



Sp 



Sp 



— —avr 



P 



I. _ dA - A— 

dv dvt dVx 



+ (^Px 

.p 



- (1 + ap 



^ _ n ^ n ^ 

dp dpt dp^ 



l-Dt^- A^ 
dp dpt dpx 



(11.60) 



—av. 



d_ 
dvt 



+ crpa 



P 



d 



d 



-{1 + apx) 
dpt dpt 



+ A 



_d_ 

dVr 



P 



+ CFpx 

P 



_d_ 

dpa 



(l + crp^) — 
dpt 



Making the changes in two last lines of Eq. (11.60) such as 



Dt 



d_ 
dvt 



d d d 

-av^Dt- Dt{avx)-^ = -av^Dt— 

dvt dvt dvt 



-Vx + avtx ] — 
P 



d_ 
dvt 



one can see that the coefficients of the differentiations (11.58) in (11.59) and 
(11.60) coincide. Inspection of the coefficients of the differentiations in higher 
derivatives Vtt, Vtx, Vxx, ■ ■ ■ requires the higher-order prolongations of the operator 
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Exercise 11.1. Verify that the operator identity (7.1) is satisfied in the same 
approximation as in Example 11.1 for the nonlocal symmetry operator from 
(11.20). 

12 Comparison with the "direct method" 
12.1 General discussion 

Theorem 8.2 allows to construct conservation laws for equations with known 
symmetries simply by substituting in Eqs. (8.23) the expressions and £ 
given by Eqs. (8.21) and (8.21), respectively. 

The "direct method" means the determination of the conserved vectors (8.3) 
by solving Eq. (8.2) for C*. Upon restricting the highest order of derivatives 
of u involved in C*, Eq. (8.2) splits into several equations. If one can solve 
the resulting system, one obtains the desired conserved vectors. Existence of 
symmetries is not required. 

To the best of my knowledge, the direct method was used for the first time 
in 1798 by Laplace [33]. He applied the method to Kepler's problem in celestial 
mechanics and found a new vector- valued conserved quantity (see [33], Book II, 
Chap. Ill, Eqs. (P)) known as Laplace's vector. 

The application of the direct method to the gasdynamic equations (11.1) 
allowed to demonstrate in [34] that all conservation laws involving only the in- 
dependent and dependent variables (11-2), (11-3) were provided by the classical 
conservation laws (mass, energy, momentum, angular momentum and center-of- 
mass) given in Section 11.1 and the following two special conservation laws 

^ J {t{p\vf + np) — px-v^dw — — Jp {2tv — x) ■ udS, 
n(t) 5(t) 

^ J {t^{p\v\'^ + np) - px ■ {2tv -x}duj^ - J 2tp {tv - x) ■ vdS 

n(t) 5(t) 

that were found in [35] in the case 7 = (n -|- 2)/n by using the symmetry ideas. 

All local conservation laws for the heat equation Ut — u^^ = have been found 
by the direct method in [36] (see in [14], vol. 1, Section 10.1; see also [37]). 
Namely it has been shown by considering the conservation equations of the form 

, X, U, tlx, "^xxj ■ ■ ■ 

, X, U, tlx, "^xx, ■ ■ ■ 

that all such conservation laws are given by 

Dt[ip{t, x)u] + Dx[uipx{t, x) - Ux(p{t, x)] = 0, 
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where v — (p{t, x) is an arbitrary solution of the adjoint equation Vf + v^x = to 
the heat equation. Similar result can be obtained by applying Theorem 8.2 for 
any linear equation, e.g. for the heat equation Ut — Au = with any number of 
spatial variables x = {x^, . . . , x"). Namely, applying formula (8.23) to the scaling 
symmetry X — ud/du we obtain the conservation law 

Dt[ip{t, x)u] + V • [u'^ifit, x) - (p{t, x)Wu] = 0, 

where v = ip{t, x) is an arbitrary solution of the adjoint equation Vt+Av = to the 
heat equation. This conservation law embraces the conservation laws associated 
with all other symmetries of the heat equation. 

Various mathematical models for describing the geological process of segrega- 
tion and migration of large volumes of molten rock were proposed in the geophys- 
ical literature (see the papers [38], [39], [40], [41], [42] and the references therein). 
One of them is known as the generalized magma equation and has the form 

Ut + [k" - u^'D^ {u-"'ut)] =0, n, m = const. (12.1) 

It is accepted as a reasonable mathematical model for describing melt migration 
through the Earth's mantle. Several conservation laws for this model have been 
calculated by the direct method in [39], [41] and interpreted from symmetry point 
of view in [42]. It is shown in [43] that Eq. (12.1) is quasi self-adjoint with the 
substitution (1.34) given by v = y}-"^-"^ iim-\-n ^ 1 and v = In if m H- n = 1. 
These substitutions show that Eq. (12.1) is strictly self-adjoint (Definition 1.2) 
if m -|- n = 0. Using the quasi self-adjointness, the conservation laws are easily 
computed in [43]. 

Some simplification of the direct method was suggested in [29]. Namely, one 
writes the conservation equation in the form (8.7), 

A (CO = . . . ,^(,)), (8.7) 

and first finds the undetermined coefficients fi" by satisfying the integrability 
condition of Eqs. (8.7), i.e. by solving the equations (see Proposition 7.1 in 
Section 7.2) 

(5 r ~ 

Ai^(a;,M,M(i), . . .) F^(a;,M,M(i), . . . =0, a^l,...,m. (12.2) 



5u 

Then, for each solution of Eqs. (12.2), the components C* of the corresponding 
conserved vector arc computed from Eq. (8.7). In simple situations C* can be 
detected merely by looking at the right-had side of Eq. (8.7), see further Example 
12.1. 
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Remark 12.1. Note that Eq. (12.2) should be satisfied on the solutions of Eqs. 
(8.1). Then the left-hand side of (12.2) can be written as 



F*{x,u,v,..., 




V=IJ.{X,U,U(^1),...) 



with F* defined by Eq. (3.2). 

The reader can find a detailed discussion of the direct method in the recent 
book [13]. I will compare two methods by considering few examples and exercises. 

12.2 Examples and exercises 

Example 12.1. (See [13], Sec. 1.3). Let us consider the KdV equation (3.6), 



and write the condition (12.2) for n = n{t,x,u). We have: 

— [fl{t, X, U){ut - U^xx - UUx)] 

ou 

d Li 

= -Dt{^J) + Dl{n) + Dx{un) - iiUx + {ut - Uxxx - ^^^)^ 

= -A(/^) + Dl{ij) + uDx{n) + {ut - Uxxx - ^'"^)^ ■ 

In accordance with Remark 12.1, we consider this expression on the solutions of 
the KdV equation and see that Eq. (12.2) coincides with the adjoint equation 
(3.7) to (3.6): 

D,{ii)^Dl{ii)+uDx{ii). (12.3) 
Its solution is given in Example 3.1 and has the form (3.8), 

II = Ai + A2U + ^3(2; + tu), Ai, A2, A3 = const. 

Thus, we have the following three linearly independent solutions of Eq. (12.3): 

/^l = 1, fi2 = U, lli = {x + tu). 

and the corresponding three equations (8.7): 



Ut — u. 



'XXX 



+ UU. 



(3.6) 



A(C^) + -Dx(C^) ^Ut- Uxxx - UUx 
Dt{C'^) + DxiC^) = u{Ut - Uxxx - u 

Dt{C') + DxiC) = (x + tu){ut - u, 




'XXX 



(12.4) 
(12.5) 
(12.6) 
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In this simple example the components C^, of the conserved vector can be 
easily seen from the right-hand sides of Eqs. (12.4)-(12.6). In the case of (12.4), 
(12.5) it is obvious. Therefore let us consider the right-hand side of Eq. (12.6). 
We see that 



(x + tu)ut — Dt [ xu + - tv? 



1 2 

2^' 



— {x-\- tu)uu^ 

-{X + tu)U:cxo. 



1 2 1 s\ 1 2 



^xx ~l~ t^x^xxj 

f 1 2 

-Dx I + 2 ^^x - ^^xx - tuu, 



Hence, the right-hand side of Eq. (12.6) can be written in the divergence form: 

{x + tu){Ut - -Uxxx - uu^) 



u 



Dt\t— + xu\+D 



Ux + t 



2 XX ^ 



— X 



u 



The expressions under Dt{- ■ ■ ) and Dx{- ■ ■ ) give and C^, respectively, in (12.6). 
Note that the corresponding conservation law 



u 



Dt\t \-xu) +D. 



u 
"3 



— X 



u 



0. (12.7) 



was derived from the Galilean invariance of the KdV equation (see [10], Section 
22.5) and by the direct method (see [13], Section 1.3.5). 

The similar treatment of the right-hand sides of the equations (12.4) and 
(12.5) leads to Eq. (3.6) and to the conservation law 



"xx g ^ ) — ^' 



;i2.8) 



respectively. Theorem 8.2 associates the conservation law (12.8) with the scaling 
symmetry of the KdV equation. 

Exercise 12.1. Apply the direct method to the short pulse equation (10.11) 
using the differential substitution (10.15). In this case Eq. (8.7) is written 



A(C^) + -Dx(C^) = UtUxt - - U^UxUxt 



U-\- - U^Uxx + UUt. \ Ut + - U^Ux -\- - U^UxUxx + 7: U^U^- 



(12.9) 
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Exercise 12.2. Consider the Boussinesq equations (8.11). Taking its formal 
Lagrangian 

L^io \M)t - gpx - fvz - iix^i^z + i^z^i^x] 
+ k + f'4^z - i>xVz + i>zVx] + r[pt + (N'^/g) i>x - i>xPz + i>zPx] , 

where uj, p, r are new dependent variables, we obtain the adjoint system to Eqs. 
(8.11): 

It is shown in [44] that the system (8.11) is self-adjoint. Namely, the substitution 

= P^-v, r^-(gyN^)p (12.11) 

maps the adjoint system (12.10) into the system (8.11). Using the self-adjointness, 
nontrivial conservation laws were constructed via Theorem 8.2. Apply the direct 
method to the system (8.11). Note that knowledge of the substitution (12.11) 
gives the following equation Eq. (8.7): 

D,{C') + D,{C') + Dz{C') 

= -0 [Axx + Azz - gpx - fVz - i^x{i^zxx + Ipzzz) + 1pz{^xxx + Ipxzz)] (12.12) 

^2 r -1 

Example 12.2. Let us consider the conservation equation (11.46), 

A(C') + D,{C^) = a{pt + vpx + pvx), 

where a is connected with the velocity v and the pressure p of the Chaplygin gas 
by Eqs. (11.22), 

1 

(7a; = , at + VGx = 0. 

P 

In this example Eqs. (12.2) are not satisfied. Indeed, we have 

/dx 
7^ 0, 
P 

[a{pt + vpx + pVx)] = o-j - Dx{(Tv) + avx = -((Jj + vax) = 0, 
[(7(pt + vpx + pVx)] = 0. 



6 

5v 
S_ 

Tp 

5_ 
5p 
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Example 12.3. Let us consider the conservation equation (11.50), 
Here Eqs. (12.2) are not satisfied. Namely, writing 

tp (^t + VV^ + ^Px^ + ^^(Pt + "^Px + P^^o;) = ^P^^t + '2.tpVVx + + tv Pt + tV^Pa; 

we obtain: 

^ [ipvt + 2tpvvx + ipa; + it'Pt + tv^px] = -t", 
— [tpvt + 2tpt't'^ + tp^ + ttipt + iv^p^] = 0. 

Exercise 12.3. Check if Eqs. (12.2) are satisfied for the conservation equations 
(11.52) and (11.54). 
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Part 3 

Utilization of conservation laws 
for constructing solutions of PDEs 

13 General discussion of the method 

As mentioned in Section 7.4, one can integrate or reduce the order of linear or- 
dinary differential equations by rewriting them in a conservation form (7.39). 
Likewise one can integrate or reduce the order of a nonlinear ordinary differen- 
tial equation as well as a system of ordinary differential equations using their 
conservation laws. Namely, a conservation law 

=0 (13.1) 

for a nonlinear ordinary differential equation 

F(x,y,y',...,y(^)) = (13.2) 

yields the first integral 

V'(x,y,y',...,yM) = Ci. (13.3) 

We will discuss now an extension of this idea to partial differential equations. 
Namely, we will apply conservation laws for constructing particular exact solu- 
tions of systems of partial differential equations. Detailed calculations are given 
in examples considered in the next sections. 

Let us assume that the system (8.1), 

Fa{x,u,U{i),...,U(^s)) ^0, Q; = l,...,m, (13.4) 
has a conservation law (8.2), 

[A(<:?0](i3.4) = O' (13-5) 

with a known conserved vector 

C = {C\ ... , C") , (13.6) 

where 

(7* = (x, u, . . .) , i — 1, . . . ,n. 
We write the conservation equation (13.5) in the form (8.7), 

A(CO = /^""F^x, u, . . . , . (13.7) 
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For a given conserved vector (13.6) the coefficients /z'* in Eq. (13.7) are known 

functions = -u, ■U(i), . . .). 

We will construct particular solutions of the system (13.4) by requiring that on 
these solutions the vector (13.6) reduces to the following trivial conserved vector : 

C = {C\x', . . . , x"), . . . , C-{x\ X--')) . (13.8) 

In other words, we look for particular solutions of the system (13.4) by adding to 
Eqs. (13.4) the differential constraints 

(x, u, . . .) = /i^(x^ x^ . . . , x"), 
(x, u, . . .) = h^{x^, x^ . . . , x"), 

C" (x,ii,ii(i),...)=/i'^(x\...,x"-i), 

where C" (x, m, . . .) are the components of the known conserved vector (13.6). 
Due to the constraints (13.9), the left-hand side of Eq. (13.7) vanishes identically. 
Hence the number of equations in the system (13.4) will be reduced by one. 
The differential constraints (13.9) can be equivalently written as follows: 

D,[C'{x,u,u^,^,..)] =0, 

Remcirk 13.1. The overdetermined system of m + n equations (13.4), (13.10) 
reduces to m + n — 1 equations due to the conservation law (13.5). 

14 Application to the Chaplygin gas 

14.1 Detailed discussion of one case 

Let us apply the method to the Chaplygin gas equations (11.15), 

P 

Pt + vp:, + pv^^O, (14.1) 
Pt + VPx - PVx = 0. 

We will construct a particular solution of the system (14.1) using the simplest 
conservation law (11.30), 

Dtip) + D^{pv) = pt + vp^ + pv^. (14.2) 
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The conservation equation (14.2) is written in the form (13.7) with the following 
conserved vector (13.6): 

C-i = p, C^ = pv. (14.3) 
The differential constraints (13.9) are written as follows: 

p = g{x), pv = h{t). (14.4) 

Thus we look for solutions of the form 

p = g{x), v^^^- (14.5) 
9{x) 

The functions (14.5) solve the second equation in (14.1) because the conservation 
law (14.2) coincides with the second equation (14.1) (see Remark 13.1). Therefore 
it remains to substitute (14.5) in the first and third equations of the system (14.1). 
The result of this substitution can be solved for the derivatives of p : 

Px = -h ^ 5- , 

9 

(14.6) 

hg' hh' h^g' 
Pt = 5-- 

g^ g g-^ 

The compatibility condition p^t — Ptx of the system (14.6) gives the equation 

g"-2^)p = g^^- 2g'h' -h^^ + 2h''C- (14.7) 

9 J h g g^ 

For illustration purposes I will simplify further calculations by considering the 
particular case when the coefficient for p in Eq. (14.6) vanishes: 

/ - 2 y = 0. (14.8) 

The solution of Eq. (14.8) is 

9{^) = — ^~~r > ciih — const. (14.9) 

ax + 

Substituting (14.9) in Eq. (14.7) we obtain 

h" + 2ahh' = 0, (14.10) 

whence 

h{t) = ktim{c- akt) (14.11) 
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if a 7^ 0, and 

h{t)^At + B (14.12) 

if a = 0. 

If the constant a in (14.9) does not vanish, we substitute (14.9) and (14.11) 
in Eqs. (14.6), integrate them and obtain 

p = k'^{ax + b) + Q cos(c — akt), Q = const. (14.13) 

In the case a — the similar calculations yield 

p^-Ax + ^AH'^ + ABbt + Q, (5 = const. (14.14) 

Thus, using the conservation law (14.2) we have arrived at the solutions 

1 



ax + 6 ' 



V — k{ax + b) tan{c — akt) , (14.15) 
p — k'^{ax + b) + Q cos(c — akt) 

and 

1 

v^b{At + B), (14.16) 
p = -Ax + ^ AH^ + ABbt + Q. 

14.2 Differential constraints provided by other conserved 
vectors 

The conservation laws (11.31)-(11.33) give the following differential constraints 
(13.9): 

pv^-p = g{x), pv + pv^ = h{t), (14.17) 
pv = g{x), p + pv'^ = h{t), (14.18) 

tpv — xp — g{x), tp + tpv^ — xpv — h(t). (14.19) 

The nonlocal conserved vectors (11.49), (11.51) and (11.53) lead to the fol- 
lowing differential constraints (13.9): 

tpv + T^ g{x), p + pv^^h{t), (14.20) 

tp^g{x), tpv-T^h{t), (14.21) 

^--s]p^ g{x), -^pv-tT^ h{t). (14.22) 
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The constraints (14.20) are not essentially different from tfie constraints (14.18). 
It is manifest if we write them in the form (13.10). 

15 Application to nonlinear equation describing 
an irrigation system 

The method of Section 13 can be used for constructing particular solutions not 
only of a system, but of a single partial differential equations as well. 
Let us consider the nonlinear equation (6.1), 

C7(V')V'* = [K{^)^.l + [K{i;) (^, - 1)], - S{i;), (15.1) 

satisfying the nonlinear self-adjointness condition (6.3), 

S'{i/j) = aC{'ip), a = const. (15.2) 

and apply the method of Section 13 to the conserved vector (9.11), 

= ,5(V')e"*, ^ aK{'il))i}^e\ = aX(V')(V'z - l)e"*. (15.3) 

The conditions (13.9) are written: 

5(V')e'^* = f{x,z), aK{'iP)'iP,e''' = g{t,z), aK{'iP){'iP, - l)e«* = h{t,x). 

These conditions mean that the left-hand sides of the first, second and third equa- 
tion do not depend on t, x and z, respectively. Therefore they can be equivalently 
written as the following differential constraints (see Eqs. (13.10)): 

aSi^jj) + S'i^lj)A^O, 

[Ki^'m, = 0, (15.4) 

[KW (V^. - 1)], = 0. 

The constraints (15.4) reduce Eq. (15.1) to Eq. (15.2). Hence, the particular 
solutions of Eq. (15.1) provided by the conserved vector (15.3) are described by 
the system 

aC{^|J)-S'{^|J)^0, 
aSi^lj) + S'i^jj)A^O, 

(15.5) 

[KWAl = 0, 
[KW (V^. - 1)], = 0. 



92 



Part 4 

Approximate self-adjointness and approximate 

conservation laws 

The methods developed in this paper can be extended to differential equations 
with a small parameter in order to construct approximate conservation laws us- 
ing approximate symmetries. I will illustrate this possibility by examples. The 
reader interested in approximate symmetries can find enough material in [14], 
vol. 3, Chapters 2 and 9. A brief introduction to the subject can be found also 
in [45]. 

16 The van der Pol equation 

The van der Pol equation has the form 

F = y" + y + e{y'^ -y') = 0, £ = const. ^ 0. (16.1) 

16.1 Approximately adjoint equation 

We have: 

^ {z [y" + y + e {y'' - y')]} = z" + z + sD, {z - 3zy'') . 

Thus, the adjoint equation to the van der Pol equation is 

F* = z" + z + e {z' - 3z'y'^ - Gzy'y") = 0. 

We eliminate here y" by using Eq. (16.1), consider £ as a small parameter and 
write F* in the first order of precision with respect to s. In other words, we write 

y" ~ -y- (16.2) 

Then we obtain the following approximately adjoint equation to Eq. (16.1): 

F* = z" + z + e (/ - 3z'y'^ + Qzyy') = 0. (16.3) 

16.2 Approximate self-adjointness 

Let us investigate Eq. (16.1) for approximate self-adjointness. Specifically, I 
will call Eq. (16.1) approximately self-adjoint if there exists a non-trivial (not 
vanishing identically) approximate substitution 

z ^ f{x,y,y') + eg{x,y,y') (16.4) 
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such that F given by Eq. (16.1) and F* defined by Eq. (16.3) approximately sat- 
isfy the condition (3.5) of nonhnear self-adjointness. In other words, the following 
equation is satisfied in the first-order of precision in e : 

= (16.5) 

Note, that the unperturbed equation y" -|- y = is nonlinearly self-adjoint. 
Namely it coincides with the adjoint equation z" -\- z — Q upon the substitution 

z = ay -\- 13 cos X + smx, a, /3, 7 = const. (16.6) 

Therefore we will consider the substitution (16.4) of the following restricted form: 

f{x,y)+eg{x,y,y'). (16.7) 



In differentiating g{x,y,y') we will use Eq. (16.2) because we make out calcula- 
tions in the first order of precision with respect to e. Then we obtain: 

z' = D^{f) + eD^{g) \y„^_y = fx + v'fy + ei9x + v'Qy - VQy'), 

z" = Dlif) + sDl^g) \^„_^ ^ + 2y'f,y + y'^fyy + y" fy (16.8) 

+ £(^^^ + 2yVx2/ - '^ygxy' + y''^gyy - '^yy'gyy' + y^Qy'y' - ygy - y'gy')- 

Substituting (16.8) in (16.3) and solving Eq. (16.5) with £ = we see that / is 
given by Eq. (16.6). Then X — C and the terms with e in Eq. (16.5) give the 
following second-order linear partial differential equation for g{x,y,y') : 

g + Dlig) I = a {4y'' - Gy'y' - 2y') 

(16.9) 

-|- (3 (sin X — 3y'^ sin x — 6yy' cos x) + 7 (3y'^ cos x — cos x — 6yy' sin x) . 

The standard existence theorem guarantees that Eq. (16.9) has a solution. It 
is manifest that the solution does not vanish because g = does not satisfy Eq. 
(16.9). We conclude that the van der Pol equation (16.1) with a small parameter s 
is approximately self-adjoint. The substitution (16.7) satisfying the approximate 
self-adjointness condition (16.5) has the form 

z ^ ay + f3cosx + 7sinx -|- eg{x, y, y'), (16.10) 



where a, (3, 7 are arbitrary constants and g{x,y,y') solves Eq. (16.9). 
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16.3 Exact and approximate symmetries 

If e is treated as an arbitrary constant, Eq. (16.1) has only one point symmetry, 
namely the one-parameter group of translations of the independent variable x. 
We will write the generator Xi — d/dx oi this group in the form (7.14): 

X. =,'!-. (16.11) 

If £ is a small parameter, then Eq. (16.1) has, along with the exact symmetry 
(16.11), the following 7 approximate symmetries ([14], vol. 3, Section 9.1.3.3): 

X2 = {^y - e [y'y' + 3xy {y' + y")]} , 

d 

X^ — {% cos X -\- e [(4 — 3y'^ — 9y^) x cos x + 3(xy^)' sin a^] } ^ , 

d 

X4 — {8 sin X + £ [(4 — 3y'^ — 9y^) x sin x — 3(xy^)' cos x] } — , 

X5 = {24y^ cosx - 2Ayy' sinx + e[[l2yy' + 9yy'^ + 9y^y')xsmx 
+ {I2y^ - 9y'^y'^ - 6y^) sin x - {l2y^ - 9y^y'^ - 9y^)xcosx 
d 

-3y^y' cos a;]}—, 

Xq^ {2Ay^smx + 2Ayy'cosx-e[{l2yy' + 9yy'^ + 9yy)xcosx (16-12) 
+ (I2y^ - 9yV^ - 6y^) cosx + (l2y^ + 9y^y'^ + 9y^)a;sina; 
d 

+ 3y^y' sin ^]}^' 
X7 = [Ay cos 2x — Ay' sin 2x + e \fi{yy''^ — y^^xcos 2x 

d 

— Sy^y' cos 2x + Gy^y'x sin 2x + 2(y — y^) sin 2x1 } — , 

dy 

X^ = {4y sin 2x + Ay' cos 2x — e [3(7/^ — yy'^^xsin 2x 

d 

+ 3y^y' sin 2x + Gy^y'x cos 2x + 2{y — y^) cos 2x\ } — • 

16.4 Approximate conservation laws 

We can construct now approximate conserved quantities for the van der Pol equa- 
tion using the formula (8.23) and the approximate substitution (16.10). Inserting 
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in (8.23) the formal Lagrangian 

C^z[y" + y + e{y"-y')] 

we obtain 

C = W[-z +e{^y''^z- z)]+W'z. (16.13) 

Let us calculate the conserved quantity (16.13) for the operator Xi given by 
Eq. (16.11). In this case W = y', W = y", and therefore (16.13) has the form 

C=-y'z' + e{3y''-y')z + y"z. 

Wc eliminate here y" via Eq. (16.1), use the approximate substitution (16.10) 
and and obtain (in the first order of precision with respect to e) the following 
approximate conserved quantity: 

C — — a (y^ + y'^) + (5 {y' sin x — y cos x) — j {y' cos x + ysm x) 

(16.14) 

+ e[2ayy'^ + 2f3y'^cosx + 2^y'^smx - yg - y'D^{g)\^„^_^j . 
Differentiating it and using the equations (16.1) and (16.2) we obtain 

Dx{C) = ey' a (%''^ — 6y^y' — 2y'^ + (3 (sin a; — ?>y'^smx — 6yy' cosx) 

+ 7 {3y'^ cos X- cos x-Qyy' sin x) - 9 - Dl{g)\^„^_^ 

where o{e) denotes the higher-order terms in e. The equations (16.9) and (16.15) 
show that the quantity (16.14) satisfies the approximate conservation law 

^-(<^)|(i6.i)^0- (16-16) 
Let us consider the operator X2 from (16.12). In this case we have 

W = Ay-e[yy + 3xy {y^ + y'^)] , 

(16.17) 

W ^Ay'-e [2y' + byy'^ + 3a; {y^ + y"")] . 
Proceeding as above we obtain the following approximate conserved quantity: 

C = 4:y\(3 cos X + J sin x) — 4y (7 cos x — (3 sin x) 

+ e[2ay'' {Ay'"" - - 2) + Ay'g - AyD,{g)\^„_^ 

(16.18) 

+ \jyy''^ — 3xy'{y^ + y'^) — 2y^ — 4y] (/3 cos x + 7 sin x) 

+ {y^y' + ^xy{y'^ + 2/'^)] (7 cos x — (3 sin x) | . 



(16.15) 
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The calculation shows that the quantity (16.18) satisfies the approximate conser- 
vation law (16.16) in the following form: 

D^{C) = 4(^ cos X + 7 sin x) [y" + y + e {y'^ - y')] 

+ Aey a (Ay'^ — 6y^y' — 2y') + /3 (sin x — 3y'^ sin x — 6yy' cos x) 

+ 7 (Sy'^cosx - cosx - Gyy'sinx) -g- Dl{g)\^„^_y + o{e). (16.19) 

Continuing this procedure, one can construct approximate conservation laws 
for the remaining approximate symmetries (16.12). 

17 Perturbed KdV equation 

Let us consider again the KdV equation (3.6), 

Ut = U^xx + UUx, (3.6) 

and the following perturbed equation : 

F = Ut — Uxxx ~ — su — 0. (17.1) 
We will follow the procedure described in Section 16. 

17.1 Approximately adjoint equation 

Let us write the formal Lagrangian for Eq. (17.1) in the form 

C^v [-Ut + Uxxx + UUx + £u\ . (17.2) 

Then 

— ^Vt- Vxxx - Dx[uv) + VUx +£V^Vt- Vxxx " UVx + SV. 

6u 

Hence, the approximately adjoint equation to Eq. (17.1) has the form 

F* = Vf — Vxxx ~ UVx + SV — 0. (17.3) 

17.2 Approximate self-adjointness 

As mentioned in Section 3.1, Example 3.1, the KdV equation (3.6) is nonlinear ly 
self-adjoint with the substitution (3.8), 

V = Ai + A2U + As{x + tu). (3.8) 
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Therefore in the case of the perturbed equation (17.1) we look for the substitution 

V = (p{t, X, u) + sip{t, X, u), 

satisfying the nonhnear self-adjointness condition 

F*| =XF (17.4) 

in the first-order of precision in £, in the following form: 

f = A + A2U + ^13(0; + tu) + eil){i, X, u). (17.5) 

When we substitute the expression (17.5) in the definition (17.3) of F*, the terms 
without e in Eq. (17.4) disappear by construction of the substitution (3.8) and 
give A = ^2 + ^3^. Then we write Eq. (17.4), rearranging the terms, in the form 

E1pu[Ut Uxxx UUx] 'i£Uxx[Ux'4'uu + '4'xu\ ^Ux[Ux'4'uuu + '^^x'lpxuu + '^'4'xxu\ 

+ e[iPt - il^xxx - uil^x + A + A2U + A^{x + tu)] = -s{A2 + ^3^)^. (17.6) 

In view Eq. (17.1), the first term in the first line of Eq. (17.6) is written e^w^u- 
Hence, this term vanishes in our approximation. The terms with Uxx in the first 
line of Eq. (17.6) yield 

ll^uu = 0, ll^xu = 0, 

whence 

f{t)u + g{t,x). 

The third bracket in the first line of Eq. (17.6) vanishes, and Eq. (17.6) becomes 

[f{t) - gx{t, x)]u + gt{t, x) - gxxxit, x) + 2[A + A^tju + A^ + A^x = 0. 
After rather simple calculations we solve this equation and obtain 

g{t, x)^Ai- Alt + (A5 + 2A2 - A^t)x, fit) ^Ae + A^t-^ A^t''. 

We conclude that the perturbed KdV equation (17.1) is approximately self- 
adjoint. The approximate substitution (17.5) has the following form: 

V Ai + A2U + As(x + tu) (17.7) 
3 



+ e 



Ae + At - ^ Ast^ )u + A4-Ait+ {A5 + 2A2 - A3t)x 
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17.3 Approximate symmetries 



Recall that the Lie algebra of point symmetries of the KdV equation (3.6) is 
spanned by the following operators: 



X - ^ X - ^ X - f^ ^ 
~ a7 ' — ^ , A3 — t- — , 

at ox ox au 

^ d d ^ d 
X4 - 3t- + x—- 2u— ■ 
at ox au 



:i7.8) 



Following the method for calculating approximate symmetries and using the 
terminology presented in [14], vol. 3, Chapter 2, we can prove that all symmetries 
(17.8) are stable. Namely the perturbed equation (17.1) inherits the symmetries 
(17.8) of the KdV equation in the form of the following approximate symmetries: 



_ d _ d 



d d d 

X4 = 3t— + X- 2u— 

at ox au 



t— - — ^s(-t^— -t— 
dx du V 2 dx du 



9 .d 



d 



t^— + 2>tx- i<otu + 2>x) 

^ dx 



2 dt 



d_ 
du 



[17.9) 



17.4 Approximate conservation laws 

We can construct now the approximate conservation laws 



(17.1) 







;i7.10) 



for the perturbed KdV equation (17.1) using its approximate symmetries (17.9), 
the general formula (8.23) and the approximate substitution (17.7). Inserting in 
(8.23) the formal Lagrangian (17.2) we obtain 



= -Wv, 

C^W [uv + V,,] - v,D,{W) + vDliW). 



:i7.ii) 



I will calculate here the conserved vector (17.11) for the operator X4 from 
1^17.9). In this case we have 



W — —2u — 3tut — xux + e ( 6tu + 3a; + - t'^ut + StxUx 1 . 



(17.12) 



We further simplify the calculations by taking the particular substitution (17.7) 
with ^2 = 1, Ai = ^3 = • • • = ^6 = 0. Then 



V = u + 2ex. 



(17.13) 
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Substituting (17.12), (17.13) in the first component of the vector (16.14) and then 
eliminating Uf via Eq. (17.1) we obtain: 

~ {2u + 3tut + xux) {u + 2sx) — e (^tu -\-2>x-\- - t^Ut + 2>txu^ u 

— + 'ituUxxx + ^tu^Ux + xuUx + eixu + dixu^xx + '^txuUx 



+ 2x'^Ux — StU^ — - t^UUxxx — t^U^U^ . 



Upon singling out the total derivatives in it is written: 



^\u^ - "ieixu + ^ tv? \ -^Dx): xu^ + tu^ - ^ tul + 3tuuxx (17.14) 

Zi \ Zi ' \- Zi Zi 

+ E {^X^U + - tXli^ — - t^U^ — 6tUx + GtxUxx + ^ t'^ul — - t^UUxx^ ■ 

Then we substitute (17.12), (17.13) in the second component of the vector 
(16.14), transfer the term Dx{. . .) from to C^, multiply the resulting vector 
(C^, C^) by 2/3 and arrive at the following vector: 



1 9 . 3 



3 2 
xu + -tu 



(17.15) 



2 

C'^ = ul — -U^ — 2uUxx + £ [xU^ — 2Ux + 2xUxx + '^tU^ — 3tti^ + QtUUxx[ ■ 

The approximate conservation law (17.10) for the vector (17.15) is satisfied in the 
following form: 

[A(C^) + DxiC^)] = 2u{ut - Uxxx - uux - eu) 

- 2£(x + 2,tu){ut - Uxxx - UUx) + o{£). 
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